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ABSTRACT 

A personalized system of Instruction combining the 
principle of auto-tutorial instruction with modified self- 
pacing was applied to a course in the fundamentals of gas 
dynamics taught to thirteen students in the Department of 
Aeronautics at the Naval Postgraduate School, Monterey, 
during the Spring quarter of 1972. The course results are 
summarized and the instructional materials developed for the 
course are included as appendices. While some degradation 
was apparent due to modifications made to the method, the 
results tended to affirm the substantial value of this 
Instructional method and its superiority to the conventional 
lecture method. 
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I. THE EDUCATION PROBLEM 



Faced with the spectre of the "post-industrial age" 

[Ref. 1] v;ith its exponential rise in population and complexity 
of life, the engineering profession is being tasked to solve 
problems of vastly increased scope v;ith reduced assets. 

Rising costs, increasingly critical public scrutiny, and a 
disturbing downv/ard trend in the number of candidates for 
engineering training are indicative of an impending crisis in 
technical education. Some of the more astute observers, 
including Vice Admiral Hyman Rlckover, have been forecasting 
such problems for more than a decade [Ref. 2]. Recognizing 
these trends, engineering educators are urgently exploring 
new methods of instruction which will not only provide more 
depth to the engineering student's education, but will also 
provide for more efficient time utilization by both students 
and teachers. 

R. Buckminster Fuller maintains that it will not suffice 
to merely expand our present educational systems, but rather 
that a "transformation" of great magnitude is essential [Ref. 
3]. The early stages of this transformation are already in 
progress and have provided some of the impetus and direction 
to this project. 

Innovations in engineering education almost invariably 
represent a departure from the venerable and currently preva- 
lent mode of teaching known as "lecturing," wherein the 
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instructor devotes a major portion of his contact hours with 
the students to the role described by Professor Plammer of 
Utah State as that of a "dispenser of Information" [Ref. ^]. 

He sees a dire need for the role of the professor to change — 
not diminishing as some observers fear with the advent of 
various other media for the dispensing of information, but 
rather assuming vastly greater Importance as the teacher 
becomes more of a diagnostician and prescrlber of learning 
experiences for each individual student [Ref. ^]. As Dr. 

Billy Koen has suggested, "... it is not unreasonable to 
expect that engineering educators ought to be able to 
•engineer' their courses." [Ref. 5] 

Learning may be thought of as resulting from the auspi- 
cious combination of several interactive factors which can 
either promote or hamper successful learning. The three 
principal Influences are motivation , operant span and stimulus 
control . Each of these is in turn influenced to a greater or 
lesser degree by several other educational factors. To draw 
an analogy from the realm of control theory, the problem in 
education is to Introduce methods which generate feedback in 
order to modify the Influence of each of these independent 
factors, so as to optimize learning. One basic premise of 
the concept of 'personalized instruction' is that this more 
effectively provides the required feedback. 

The influencing factors are defined as follov/s: 

Motivation - The force which drives a student to seek to 
learn . 
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Stimulus Control - The effect exerted by the learning 

surrounds on the Individual. 

Operant Span - The maximum scope (In size or complexity) 

of a single block of material v/hlch may be 
assimilated by a given student at a time. 

Structure - The concept and method of presentation of the 
course material. 

Facilities - The physical surrounds (Including graphic aids 
and other media). 

Pacing - Rate of delivery. 

Intellectual Ability - The student's Inherent talent. 

Level of Development - The student's specific level of 

achievement or skill In a given 
subject area. 

Quotidian Variability - A psychometric term used to Identify 

dally variations In the student, not 
Intrinsic to the educational method, 
such as Illness, fatigue, or behavior 
arising from Idiosyncratic response 
to some unique aspect or feature of 
the educational environs. [Ref. 6] 

The Interaction of these factors forms the basis of a 

substantial portion of current educational research. Detailed 

Investigation of these Interactions represents subject matter 

well beyond the scope of this paper. Structure and pacing 

are the most obvious areas of concern In a system of 

personalized Instruction , though each factor must be considered. 

In the development of better methods, two unique aspects 

of engineering education must be considered: 

a. The professional skills and concepts Involve principally 
the cognitive domain. 

b. A hierarchy of skills exists, which causes mastery of 
most areas to be dependent on a thorough grounding In 
preceding material. 
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Failure to thoroughly master the required hierarchy of 
skills (a common tendency among many students) can lead to 
the development of a condition referred to as "cumulative 
Ignorance," which probably accounts for the majority of the 
failures and "drop outs" In the engineering educational 
process. If no other gain be realized, the possibility of 
precluding some cumulative Ignorance justifies the efforts 
expended In searching for better methods. However there Is 
mounting evidence that methods have been devised and are 
proven which not only help solve the problem of cumulative 
Ignorance, but also provide far more optimal Interaction of 
the various parameters affecting the learning process. 
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II. THE CONVENTIONAL LECTURE METHOD 



The standard with which all other educational methods are 
usually compared (or contrasted) is the "conventional 'lecture 
method." The term as used herein connotes a formal discourse 
wherein the instructor ("lecturer") presents the course 
material through the media of the spoken word, using no 
materials or training aids other than a chalk board. The 
lecture method normally exhibits the following characteristics: 

a. The lecture is time-constrained by the standard 
scheduled period. 

b. The student is in a passive role and is expected to 
transcribe by notes, or otherwise record, whatever 
information he deems essential. 

c. The pace is essentially fixed by the lecturer based on 
the scheduled contact hours/course coverage. 

d. There is virtually no "real-time" feedback. Due to the 
time constraints, questions or individual difficulties 
are usually deferred. This leads to the not unusual 
situation of the student missing vital succeeding 
concepts while momentarily pondering a puzzling point. 

e. The clarity of presentation is strongly dependent on 
the expository ability of the lecturer. 

The conventional lecture method is seldom utilized exactly 

as described above. It is normally modified by one (or more) 

of many variations in structure or facilities ranging from 

simple devices, such as the use of colored chalk, to the comple 

such as closed circuit television with printed, pre-distributed 

notes. V/here class size and time permit, most lecturers 

seek feedback by encouraging questions from the class at 

any time. While these variations mitigate somev/hat the 

shortcomings of the lecture, the basic vjeaknesses remain. 



The more capable students are held back and the students 
suffering from any substantial amount of "cumulative Ig- 
norance" find their burden of cumulative Ignorance Is grovjlng. 

The lecture method has three advantages vjhlch have acted 
to foster Its acceptance for hundreds. If not thousands of 
years. On closer examination these 'so-called* advantages 
prove to contribute to administrative expediency - but often 
at the expense of learning. There Is truth In the statement 
that students often learn the material of a course In spite 
of the teaching methods rather than because of them. Firstly, 
the conventional lecture provides for the most rapid and 
Inexpensive output of Information by the lecturer with the 
minimum amount of equipment or materials. This says nothing 
about the amount of Information actually received by the stu- 
dent. Secondly, the elegance of presentation and the Influence 
of the lecturer on his class can serve to "sell" the concepts 
Involved to the class. Thirdly, this method Is accepted. 

This Is a very real consideration v;hen one must deal v;lth 
the problem of Institutional Inertia. 
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III. ALTERNATIVES TO THE CONVENTIONAL LECTURE METHOD 



Many of the variations or innovations associated with 
non-lecture methods, such as v/rltten course objectives, 
instructor tutoring, peer tutoring, visual aids and printed 
notes or outlines, can be effectively used to enhance the 
lecture. However, while improving the lecture method, these 
innovations do not directly contribute to the basic need of 
each student to explore the topic ^ his own best speed . 

This can only be achieved by abandoning the lecture as the 
primary information vehicle. 

Currently in vogue among engineering educators are 
several varieties of non-lecture methods, nearly all based 
on the reinforcement theories first expounded by Sidney 
Pressey and B. P. Skinner [Ref. 7]. Their research resulted 
first in the development of "programmed Instruction". The 
Programmed Text and the Teaching Machine are direct applica- 
tions of this reinforcement theory and are now widely used 
(vjith proven effectiveness) for the instruction of material 
involving verbal or skill knov/ledge. These methods suffer 
from relatively Inflexible structure and the tendency to 
design the Instructor out of the picture. Also, neither 
method has proven quite so adaptable to complex concepts 
involving cognitive strategy, such as engineering problem 
solving or design. For this reason, engineering educators 
have been among the last to experiment with nevi teaching 
techniques. However, these shortcomings have recently 
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prompted the development of several closely related methods 
knovm variously as "Auto-Tutorial Instruction" [Ref. 8], 
"Self-Paced Individually Prescribed Instruction" [Ref. 9], 
"Personalized (or Proctorial) System of Instruction" [Ref . 10, 
11] or some similar designation, all of v;hich are more or 
less self-paced and vary primarily in the method by which 
reinforcement is achieved and progress is evaluated. However, 
in all cases the teacher , in the role of diagnostician, tutor, 
counselor and designer of learning experiences [Ref. 12] is 
the key to success of these systems. 
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IV. PERSONALIZED INSTRUCTION 



Professor P. H. Keller is generally regarded as the 
principal architect of this method which he calls the 
"Proctorial (or personalized) System of Instruction" (PSI), 
though many of the features were concurrently and indepen- 
dently developed by Professor S. H. Postlewaite at Purdue 
University in Indiana. The basic features of this method 
are succinctly described by Professor B. V. Koen of the 
University of Texas as follovjs: [Ref. 13] 

1. The ' go-at-your-own-pace ' feature, v/hich permits a 
student to move through the course at a speed commensu- 
rate v;ith his ability and other demands upon his time; 

2. The unit-perfection requirement for advance, v/hich lets 
the student go ahead to nevj material only after 
demonstrating mastery of that vjhlch preceded; 

3. The use of lectures as vehicles of motivation, rather 
than sources of critical Information; 

The related stress upon the written work in teacher- 
student communication; and, finally: 

5. The use of proctors, which permits repeated testing, 
immediate scoring, almost unavoidable tutoring, and 
a marked enhancem.ent of the personal-social aspect of 
the educational process. 

This summary by Dr. Koen clearly Illustrates the particular 
desirability and applicability of this method to engineering 
education and accounts for the growing success of this method. 
The particularly agonizing problem of "cumulative ignorance" 
is effectively overcome by the unit-perfection requirement of 
this method. 

The concept of PSI is predicated on the validity of 
* reinforcement theory . ' The course is divided into an 
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appropriate number of units, depending on major subject areas 
and complexity of the material. Each student is provided 
with a Study Guide for the unit of instruction on which he is 
working which will contain: 

a. Course objectives - defined in terms of measurable 
behavior. 

b. Reading assignments (if a text is used). 

c. Discussion of theoretical concepts and application 
to specific situations. 

d. Example problems with a discussion of the general 
philosophy and specific techniques of solution. 

e. Practice problems, with ansv;ers and guidance vjhere 
necessary . 

The student's initial exposure to new concepts can thus proceed 
at his own pace. A fundamental concept of PSI is that the 
instructor holds the key to success of the system in his role 
as "trouble-shooter." Ideally the Study Guides will provide 
a sufficiently clear exposition of the concepts so that the 
teacher will be freed to spend most of his time dealing with 
the individual problems which require personal attention. 

The principal problems encountered by practitioners of 
this method were summarized recently by Dr. Gilmour Sherman 
of Georgetown University — vjho, as an associate of Professor 
Keller in the development of the PSI method, has received 
considerable feedback from various educators implementing PSI 
in one or another form. Dr. Sherman's list can be summarized 
as follows: [Ref. l4] 

1. Professor works harder than ever. 

2. Students work harder. 

3. Substantial increase in logistic and administrative 
workload . 



14 



A. 



. staff (if paid) becomes expensive. 

5. Class size for effective implementation is limited. 

6. Production lead time for materials is burdensome. 

7. Procrastination by students creates problems with 
the administration. 

8. Hovj to assign grades? 

Of these, the last two problem areas seem the most pernicious. 
The advocates of PSI with experience have generally shovm 
that the system can be engineered to minimize all of these 
problems with a few iterations, resulting in proven increases 
in effectiveness of Instruction. For an in-depth current 
assessment of the Keller Method of PSI, the reader is referred 
to Ref. 15, in its entirety. 
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V. APPLICATION OF PSI TO A GAS DYNAMICS COURSE AT NPS 



An Initial consideration in the development of PSI 
materials for the Gas Dynamics course offered at NPS was 
to insure compatibility with the quarter system. It was 
therefore determined that while retaining the self-paced, 
self-guided features of the course for each instructional 
unit, the pace would be controlled in a longterm sense by 
scheduling the check tests (to determine mastery of each 
unit) to be taken at regular Intervals by all ( 13 ) students 
involved in the PSI course. This also permitted adequate 
lead time for the development of the Study Guides for each 
unit and allowed some feedback of Information to modify 
later units to Improve any notable weaknesses in format 
discovered in early units. The study guides are included 
as Appendix A. 

Tv;o major mid-course modifications were instituted. The 
first, a decision to discontinue the practice of repeating 
check tests (except in an extreme case) was made by the 
instructor when it was determined that time constraints may 
have prevented the class from finishing within the quarter. 
Secondly, in order to appease a number of students who found 
it discomforting to have no group sessions, the instructor 
instituted the practice of group review discussions in order 
to re-lnforce the material. This provided tutoring on a 
group basis, for those who desired to attend the problem 
sessions . 
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VI. EVALUATION OF RESULTS 



The success of this project was Intended initially to 
be measured by three principal means: 

1. Comparison v^;ith a "control" group in the form of a 
class of tvjelve students taught by the instructor 
using the conventional lecture method during the 
same quarter. 

2. A comprehensive student quest ionaire administered 
at the end of the course to the thirteen students 
Involved in the "PSI" class. 

3. A critique session concerning the course, to be 
moderated by the writer at the end of the quarter. 

Due to the lack of statistical significance of any data 
involving groups of this size, grade comparisons are 
considered to have little meaning. However the instructor 
felt that there vjas a tendency toward better overall 
comprehension by the students in the PSI group as determined 
by major quiz and final examination averages. 

The questionalre (Appendix B) vjhlch was given to each 
student in the PSI group revealed several key facts. Tv/elve 
of the thirteen students returned the questionalre and 
where appropriate the numerical average is tabulated in the 
right hand column beside the question. Notev;orthy are the 
student's apparent high regard for the concept of defined 
objectives, study guides and check tests. Significant 
comjnents from the last questions are included in Appendix C. 

The critique session tended to corroborate the comments 
in the quest lonaires . A notevjorthy comment from this session 
is included in Appendix C. 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



In-so-far as this course followed the model of PSI as 
conceived by Professor Keller, the results, based on the 
subjective evaluation of the students involved and the 
instructor, tend to substantiate the value of this approach. 
Given the bastardization of the long-term self-paced concept 
which resulted from controlling the administration of the check 
tests in the manner described, the hoped-for uniform high 
level of mastery was not achieved. However, several comments 
by students in the course indicated the ’depth* of knowledge 
achieved was generally greater than normal and most students 
admitted they vjorked harder than in previous courses. The 
controls applied to the pacing of the course negated a true 
evaluation of the problem of student procrastination; however, 
comments from tvjo students who fared poorly in the course 
attributed this to a personal lack of interest in the subject 
area. It is also notev.^orthy that only four of the twelve 
students who returned questionaires sought tutorial aid from 
the instructor. This resulted in part from the students' 
dependence on the review discussion sessions to iron out 
problems and in at least one instance was the admitted result 
of a student’s personal antipathy tov/ard the course and 
instructor. This is a fairly obvious example of quotidian 
variability adversely affecting motivation, hence learning. 

The student in question is otherwise average but fared poorly 
in this course. 
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It Is recommended that in the next offering of Gas- 
Dynamics the course be made truly self-paced with the use of 
a proctor drawn from one of the students v-zho have taken the 
course under this project. Eleven of the twelve students 
who responded to the questionaire indicated their willingness 
to serve as a tutor. While it is not unreasonable to assume 
that any student volunteering to serve in this capacity would 
be substantially compensated by the challenge Involved in 
this task, the time constraints and the inherent demands of 
the academic load certainly justify some more tangible reward. 
Many institutions are paying the graduate students utilized 
in this capacity. This being inappropriate at this institution, 
it is strongly recommended that some sort of academic credit 
be granted to the student undertaking this task as tutor or 
proctor . 

\7ith the materials and experience resulting from this 
first cycle, the foundation exists to truly exploit this 
proven method in succeeding offerings of the Gas Dynamics 
course at the Naval Post-graduate school. Moreover the 
flexibility and effectiveness of this method indicate the 
desirability of applying this method to other courses. 

Perhaps the most concise statement to date concerning the 
underlying strength of the Keller method was made by 
Professor A. J. Dressier in the foreword to a recent work- 
shop on the Keller Method. [Ref. 15] 
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"The rapid spread of the Keller method is remarkable 
when one considers the handicaps it has had to over- 
come : 

(1) It is a new idea that challenges a teaching 
technique (the lecture method) that is hundreds, if 
not thousands, of years old. New ideas are not usually 
received with delight. 

(2) It is based on a controversial theory; the 
reinforcement theory associated with the name of 

B. P. Skinner. 

(3) It Involves more v/ork for both the teacher 
and the student. More work is something that most 
people feel they do not need. 

In spite of these negative features, the Keller method 
has spread at a remarkable rate for one simple reason - 
it works!" 
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APPENDIX A - STUDY GUIDES 



The Study Guides written for this course were developed 
to be used in conjunction with the text Gas Dynamic s by 
Professor James E. A. John [Ref. l6]. Additional reference 
material and assistance vjas gained from Refs. 1? - 20. 

The Course Objectives were developed using the methods in 
Refs. 21 and 22. 

Some additional handouts were prepared by the Instructor 
and distributed to assist the students with certain points. 
Some of this material will eventually be incorporated in 
the future revisions of the Study Guides. 
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INTRODUCTION 



It is assumed that before entering the world of Gas 
Dynamics you have had a reasonable background in Mathematics 
(through Calculus), together with an exposure to basic Fluid 
Dynamics and a course in elementary Thermodynamics. 

Prior to undertaking nev; material it will prove fruitful 
to review some fundamental concepts and facts from which 
the study of Gas Dynamics will proceed. 

First read Sections 1.1 through 1.5 in the text. You 
may also find it useful to review your notes or textbooks 
used in previous Fluid Mechanics and Thermodynamics courses. 
Then answer the following review questions as completely as 
possible using any source you wish. Your answers are to be 
submitted to the instructor at which time you viill be pro- 
vided correct answers to these review questions. No attempt 
should be made to continue with this course until you feel 
that you understand the answers to all of these review 
questions . 



* All references to sections, equations, diagrams, etc. in 
the 'Text' refer to GAS DYNAMICS by James E. A. John. Any 
other referenced sources v;ill be referred to by name. 
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REVIEVJ QUESTIONS 



1. Can you rocall tho Taylor’s Serios Expansion, It s 
applications and limitations? 

2. Kow is a derivative such as ^ defined? 

3. Cite the basic units of measurement in the MET system, 

(Metric and English) 

4 . VJhat is the definition of a one pound force (in terms 
of pounds-mass, feet and seconds)? 

5. Explain the significance of g^ in Newton's Second Law. 

6. V/hat are a consistent set of units for g^. 

7. V.’hat is the distinguishing characteristic of a fluid 
(as compared to a solid)? 

8. Describe the basic difference between gases and liquids. 
VJhen can a gas be treated as a continuum? 

9 « Explain the difference betv/een a microscopic and a 

macroscopic approach to the analysis of fluid behavior. 

10 . Describe a FLUID PROPERTY and give three or more examples 
of a fluid property. (Note: We will use the term Fluid 
property’ to mean ’Gross’ or ’Bulk’ fluid property as 
used in the text . ) 

11 . Fluid properties may be categorized as either INTENSIVE 
or EXTENSIVE. Define what is meant by each and list as 
many examples of each type of property as you can. 

12 . V.Tiat does the term "specific” imply? (example: specific 
enthalpy ) 

13. What is the relationship betv^een density and specific 
volume? 

14 . Describe an EQUATION OF STATE and recall one with which 
you are familiar. 

15. State the PERFECT GAS LAW. When is it valid? (i.e. what 
are it’s limitations?) 

16. Give a consistent set of units for each term in the 
Perfect Gas Law. 

17. Define a system. Describe the difference between an 
open system and a closed system. 
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18 . Describe the CONTROL VOLUME APPROACH to a problem and 
compare It with the CONTROL MASS APPROACH. Can you 
identify the INDEPENDENT and DEPENDENT variables in 
each case? 

19. Define a STATE FUNCTION and a PATH FUNCTION. Give 
examples of each. 

20 . V/hat is a PROCESS? ' What is a CYCLE? 

21 . How is the Zero th LAW of THERI-IODYNAMICS related to 
temperature? 

22 . State the FIRST LAW of THERMODYNAMICS for a closed 
system. 

23. Give a form of the FIRST LAW applicable to an open 
system. 

2 k. State any form of the SECOND LAW of THERMODYNAMICS. 

25.. Define ENTROPY. 

26. Give equations that relate entropy to 

a. ) Internal energy. 

b. ) enthalpy. 

27. Define a REVERSIBLE PROCESS for a THERMODYNAMIC system. 
Is any real process ever reversible? Of what practical 
value is the concept of reversible processes? 

28. What are some IRREVERSIBLE EFFECTS? 

29. V/hat is an ADIABATIC process? 

30. If a process is BOTH REVERSIBLE and ADIABATIC it is 
also something else. What? 

31. Vdnat STATE VARIABLE is constant in an ISENTROPIC 
PROCESS? Is it possible to have an Isentropic process 
that is not a reversible adiabatic process? 



32. If we say dS = 0 , does this automatically imply an 
ISENTROPIC situation? (where "S" is ENTROPY). 

33. Define (in the form of a partial derivative) the SPECIFIC 
HEAT'S (Cy and Cp). Are these expressions valid for 

a material in any state? Are they valid for any process? 

3 ^. For a PERFECT GAS, SPECIFIC INTERNAL ENERGY is a 
function of which state variables? 

35. For a PERFECT GAS, SPECIFIC ENTHALPY is a function of 
v.’’hich state variables? 
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36. For the equation: Q = W + AE 

How are the signs conventions defined for Heat and V/ork? 

37. Process Plots 



P 



n increasing 




V 




pv^ = ; where (|^ is a "constant" 

n is any positive number. 

This is a General Polytropic Process. 

a) In the p-v diagram place each of the following in 

the correct box: 0,(p , S. 

b) In the T-s diagram, label each process line with 
the correct value for n and identify which fluid 
property is held constant. 

c) Are these plots general, or are they restricted 
to a Perfect Gas? Why? 
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COMMENTS OF THE REVIEW QUESTIONS 



1. Prom CRC Tables (17th edition) 

TAYLOR'S SERIES; 

2 

f(x)=f(a)+ (x - a)f'(a) + ■ f"(a) + 

The Taylor's Series expansion enables one to examine the 
unknown fluid properties at an arbitrary point (x) in a 
continuous flow field by an expansion of the derivatives 
of the properties at a point (a) where these are known. 
Note that all of the derivatives are evaluated at point 
(a). 



dy _ lim f(x + Ax) - f(x) 
dx Ax->0 Ax 

The important point here is that the derivative is a 
process which becomes valid in the limit as the incre- 
mental change in the state variable becomes vanishingly 
small . 



3 . SYSTEM 





(MKS) 


Absolute 

Metric(CGS; 


British 
Gravitat' nl 


Absolute 

English 


English 

Engineering 


Mass 


kilogram 


gram 


slug 


pound- 

mass 


pound- 

mass 


Length 


meter 


centimeter 


foot 


foot 


foot 


Time 


second 


second 


second 


second 


second 


Force 


newton 


dyne 


pound- 

force 


poundal 


pound- 

force 


Sc 


kg-m 

nevjt on-sec 


^ gm-cm 


^ slug-ft 


^ Ibnrft 


32.17^ 

ft-lbm 

Ibf-sec^ 


^ dyne-sec-^ 


lbf-sec2 


poundal- 

sec^ 



A one pound force will give a one pound mass an accelera- 
tion of 32,17^ ft/sec2. 

5. FORCE = MASS times ACCELERATION 

FORCE = (Constant ) (Mass ) (Acceleration ) 
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The CONSTANT rectifies the dimensions Involved. If we 
assign the symbol "1/g^" to this constant, Nev;ton's Lav; 
can be stated 




and g^ will have the values Indicated In the table above. 



6 . 



1 (Ibm) 32.17^ (ft/sec^) 
®c 1 (Ibf) 



32 . 17 ^ 



ft-lbm 

lbf-sec2 



7. It v;ill continuously deform under the influence of an 
applied SHEARING FORCE. 

8. There is no rigorous distinction between gases and liquids 
but generally we think of liquids as being nearly incom- 
pressible (actually they are slightly compressible) v;hereas 
gases are very compressible (thus gases assume the volume 
of a closed container). The problem in making a clear 
distinction becomes impossible when one considers a 
substance which is above its critical point. 

9. Refer to Section 1.2 in the text ( Gas Dynamics , John) 



10. Considering the macroscopic approach of a small, but finite 
volume of fluid consisting of a sufficiently large number 
of molecules, we may define a FLUID PROPERTY by using the 
limit concept of two relationships such as: 

p _ lim / normal >. where "a" is a small 

AA-^a ^ AA but finite area. 

In any case the resulting characteristic exhibited by the 
fluid will be as though we were considering a single point 
in the fluid. Such a characteristic (a property) is 
independent of the path by which the point was reached. 

Examples are: pressure, density, tem.perature , entropy, 

enthalpy, etc. 

11. Intensive properties exhibit the same characteristics 
Independent of the amount of mass present. Temperature, 
pressure and specific properties are INTENSIVE. 

Extensive properties are dependent on the amount of 
mass present. Mass, energy, entropy and enthalpy are 
examples . 

12. The term "specific" refers to the amount of an extensive 
property per unit mass. 
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13. Density Is the reciprocal of specific volume: 
1 



1^1. An equation of state is a relationship among properties. 
The most familiar is the Perfect Gas Law: 



P = pRT 



L - pi 






where P 


is 


pressure 


P 


is 


density 


R 


is 


a gas const 


T 


is 


temperature 



15. The Perfect Gas Lav; is derived from kinetic theory and 
neglects molecular volume and inter-molecular forces. 
Therefore, the Perfect Gas Law can be used to deal with 
gases under conditions of low density (consider the 
variables) v;herein these two phenomena are not significant 
factors. Gases at lov; pressures and high temperatures 
are closely described by the Perfect Gas Law. 

16. P — Ibf/ft^ (absolute pressure) 
p — lbm/ft3 

T — °R (absolute temperature) 

R _ ft-lbf 
^ lbm-°R 

17* The term system is used in a broad sense to identify 

clearly whatever it is that is under consideration. A 
closed system is one in which mass cannot enter or leave; 
an open system is one in which mass may enter and/or 
leave . 

l8. The CONTROL MASS approach to a fluid flov; problem utilizes 
a closed system, whereas the CONTROL VOLUME approach 
considers a fixed volume in space. Energy can cross the 
boundaries of either system v;hereas mass can only cross 
the boundaries of an open system. 

19- A state function (also called a point function or property) 
is an observable characteristic which describes the 
physical or thermal state of a system. Tv;o or more state 
functions fix the state of a system. 

Path functions are transient phenomena. Consider Heat 
and Work: 

— Systems never possess heat or work, but either or both 
cross the system boundary when a system undergoes a 
change of state. 
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— Both heat and v;ork are boundary phenomena. Both 
are observed only at the boundaries of the system, 
and both represent energy crossing the boundary 
of the system. 

— Both are PATH functions and inexact differentials. 

20. A process is a change of state in a system which has a 
unique sequence or relationship between the state 
variables. A cycle consists of several processes which 
vary the state of a system and return it to the original 
state . 

21. The so called "zeroth" law of thermodynamics defines the 
state variable temperature in terms of tvjo systems which 
are in thermal equilibrium with a third system. 

22. For cyclic operation: ZQ = EW or ^dQ = 

For non-cyclic operation: Q = W + AE 

6Q = 6W + dE 

where heat flow into the system is positive and 
work transferred out of the system is positive. 

Thus a new property has been defined — ►E, the total 
energy of the system. 

k / « » 'J'' + / ^ 0 (v-dA) 

C . V c . s . 



2^. It is impossible to make any transformation whose only 
final result is the exchange of a non-zero amount of 
heat vjith less than two heat reservoirs and the appear- 
ance of a positive amount of work in the surroundings, 
(from Lord Kelvin's statement of the second law of 
thermodynamics, I 85 I) 



25. 

26 . 



dS 



6Q 



REV 

T 



For any process: 



Tds = du + pdv 

where "u" is specific internal energy, "s" is specific 
entropy and "v" is specific volume. T & p are 
temperature and pressure. 

Tds = dh - vdp 

where "h" is specific enthalpy. 



% ^ 

INTRODUCTION TO THERMODYNAMICS, Sonntag & Van Wylen. 
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27 . Conceptually, a reversible process is one in vjhich both 
the system and its surroundings can be returned to their 
original conditions. Although NO real process is 
reversible in a strict sense, processes which closely 
approximate this ideal may be realized. 

28. An irreversible effect causes a disturbance in the 
equilibrium state of a system. To quote from Elements 
of Gas Dynamics , by Liepmann and Roshko; 

... A system is in equilibrium if it is free of 
"currents". The term current refers to the flux of 
a quantity like heat, mass, momentum, etc. A current 
of heat flows if there is a finite temperature 
difference; a current of mass flows if there is a 
difference in concentration of one component; a 
current of momentum flows if there is a difference 
in velocity. 

Thus stirring a fluid, sudden heating or any process 
involving friction will introduce irreversibilities. 

29. A process involving no transfer of heat to or from the 
system is called adiabatic. 

30. Isentropic. 

31. Entropy. It is possible to balance the increase of 

entropy xvithin a system due to irreversibilities with a 

decrease in entropy due to heat flow and achieve a 

resultant change in entropy of zero (or an ISENTROPIC 

process). This process is obviously neither adiabatic 

nor reversible. 

Recall; dS = dS + dS. 

e 1 



32. Yes. 

33. 

These definitions are valid for any process and for any 
material in any state. 

34. Temperature ONLY. 

35. Temperature ONLY. 

36. Heat into the system is positive. Work done by the 
system is positive. 



= ( 



9 T 
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37. 





V s 



The process lines with n = 1 referring to T = 
constant ((^) and n = y referring to isentropic are 
restricted to Perfect Gases. V/hy? 
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UNIT 1 - BASIC EQUATIONS - OBJECTIVES 



The student shall be able to: 

1. State the basic concepts from which the study of Gas 
Dynamics proceeds. 

2. Explain the difference betv/een a MATERIAL DERIVATIVE 
(or TOTAL DERIVATIVE) and a PARTIAL DERIVATIVE with 
respect to time. 

3. Obtain the properties at one point in a flow field in 
terms of properties and their derivatives at another 
point by means of a Taylor's Series Expansion. 

4. State the equation used to relate the material derivative 
of any extensive property to the properties inside of and 
crossing the boundaries of a control volume. Interpret 
each term in the equation. 

5. Starting with the basic concepts or equations which are 
valid for a CONTROL MASS, obtain the integral forms of 
the ENERGY, CONTINUITY, and MOMENTUM equations for a 
CONTROL VOLUME. 

6. Simplify the integral forms of the ENERGY, CONTINUITY 
and MOMENTUM equations for a control volume for the 
conditions of steady, one-dimensional flov/. 

7. Apply the simplified forms of the ENERGY, CONTINUITY, 
and MOMENTUM equations to differential control volumes. 

8. Show that by introducing the concept of ENTROPY and the 
definition of ENTHALPY, the path function HEAT (dQ) may 
be removed from the ENERGY equation to yield the expres- 
sion sometimes called the 'PRESSURE ENERGY Equation.' 

+ Sd2 + TdS. + 6W =0 
P 2g^ g^ 1 s 

9. Simplify the PRESSURE ENERGY Equation to obtaan BERNOULLI'S 
Equation . 

10. Apply Nevjton's Second Law of Motion to a differential 
control volume and develop the MOMENTUM equation for 
Steady, One-dimensional flow. 
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11 . 



Identify each term and explain the restrictions and 
assumptions for each of the follov/ing equations: 



1 ) 


6^ 

p 


+ 


dA 

A 


II 

0 




2) 


6 q 


= 


d(Pv) 


+ < 5 v; . 

shaft 


+ du + 


3 ) 


P 


+ 


d(y2) 

2gc 


+ gdz ^ 
Sc 


Tds^ + 




P 


+ 


f (V^) 
‘’e^So 


dx + 

Sc 


VdV 

Sc 

2 


5 ) 


£1 

P 


+ 


2gc 


II 


P2 V2 

P 2g^ 



d(V^) 

2gc 



6W3 = 

0 

+ g^2 

Sc 



+ 

0 



gdz 

Sc 
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study Guide Unit 1 
Basic Equations 



1.1 Equations of Motion - 

a) Reviev/ Section 1.^ of the Text. List the unknowns and 
the Equations (by name) v;hich v;ill be used to solve 
for these unknowns in the study of compressible fluid 
flow. How does the problem simplify if the flow is 
considered incompressible? 

b) State in your own words v;hat is meant by (1) the 
conservation of Mass, (2) the Conservation of Energy 
and (3) the Conservation of Momentum. Is Momentum 
actually conserved? 



The Basic Relationships with vjhich we will analyze Gas 
Dynamics will be: 

(1) CONTINUITY (Derived from Conservation of Mass) 

(2) ENERGY (Derived from Conservation of Energy) 

(3) MOMENTUM (Derived from NEWTON'S SECOND LAW) 

(4) An EQUATION OF STATE p = (|: 

p = pRT 
others 

(5) THERMODYNAMIC CONSIDERATIONS 



3 ^ 






I 
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1 . 2 Transformation of a Material Derivative of any Extensive 
Property to a Control Volume Approach. 



Let X be the total amount of any extensive property In 
a given mass. 

Let X be the amount of X per unit mass. 



Thus : X = J' xdm = fff p xdv — J' p xdv (*1.21) 



V 



where dv e Incremental volume element. 




Consider what happens to a material derivative such as 



dt 



dt 



11m 

At->0 



(Final value of ~ (Initial value of X)^ 

At 



(* 1 . 22 ) 



Final X X of regions 2 & 
Initial X -> X of regions 1 



3 computed at tlme^t + At* 

\ ! 

& 2 computed at time t ~ 



dt 



11m 

At-^0 




xdm + 




t + At 
At 




+ 




(*1.23) 
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Consider the two terms: 



lim 


7 


^ xdm 1 


’ 1 


xdm 


At->0 




? Jt+At 


J‘ 


2 



At 



But since 




lim 
A t"^0 



+ At) - X2(t) 
At 



8X2 

9t 



1. Partial derivative notation is used since the region is 
fixed . 

2. As At Oj region @ the original confines of the mass. 
We shall call this region the control volume. 



Thus : 



lim 

At->0 




t + At 
At 




xdm 



t 



9X 

c 



V . 



9t 




pxdx (*1.2^ 



Consider the term 



lim 

At->0' 



f xdm 




/3 J 


t + At 



At 



Note - by definition 
region (D is formed 
by the fluid moving 
out of the control 
volume . 



This integral represents the amount of X in region (J) 
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Let:n be a unit outward normal 



Note: n ds Is equivalent to dA in the Text, 

Section 1.6. 



(V.n) 
p(V.n) 
p(V.n) ds 

^ A* 

px(V.n) ds 



ds be an Increment of surface dividing 
regions (2) & (3) 

V.n = component of V to ds 

ds = volumetric flow rate 

ds = mass flow rate 

At = amount of mass that crossed ds in time At 
At = amount of X that crossed ds in time At 



Thus 



out 



px(V.n)dsAt ~ total amount of X in region 3 ( 



-1.25) 
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and 



lim 
At "^0 



t + At 



(* 1 . 26 ) 



J ' sdm 



At 




xp (V. n) 



ds 



where S , is the surface vjhere fluid leaves the control volume, 
out 



WHY THE LIMIT? All properties are going to be evaluated at 

the surface S, thus equation (*1.25) is only 
an approximation. It becomes exact in the 
limit At->0. 

This integral is called a flux or rate of X flow out of the 
control volume. 



Similarly : 



with n' a unit inward normal 




Note that this term appears in equation (*1.23) with a minus 
sign. If we use n (outward) here, since n = -n', the 
combination of these two terms can be obtained by integrating 
the following over the entire surface: 



S 



/ 

out 



xp(V.n) ds 




(V.n' ) 



ds 




xp(V.n) ds 



(*1.27) 



v;here c.s. represents the control surface surrounding the 
control volume. 
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Pinal result becomes: 



dt material 

derivative 



9t 




Triple 



xpdv + 



integral 




xp(V.n)ds (*1.28) 
Double integral 



AGAIN: This applies to any extensive (Mass dependent) property 

vjhether it is a scalar or a vector. 



In VJords: The rate of change of X for a given mass as it is 

moving around is equal to the rate of change of X 
Inside the control volume plus the net flux (flow 
rate) of X out of the control volume. (i.e., vjhat 
goes out minus what comes in). 



X can be momentum -> x = V (used to transform Newton's 2nd law) 



X can be energy x = e = u + ^ 

C 0 

(used to transform 1st law of thermo) 



X can be mass itself -> x = 1 (used to transform continuity eq.) 



or other. 



Example: This may help to sort out v/hat is and is not a 

Material Derivative (or the change of an extensive 
property within a control mass) 

For a Control Mass: 




These are NOT Material 



Derivatives as they des- 
cribe the rate of transfer 
of energy across a surface . 



/ 



dt 

This ^ a Material 
Derivative . 
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Thus we may summarize the tv;o approaches to a flow problem: 



1. ) CONTROL MASS - wherein a given mass is observed. 

This is called a Closed System . (T-e., no new 
matter can enter or leave the system) 

2. ) CONTROL VOLUME - wherein a particular fixed volume 

in space is observed. This is called an 
Open System. 



As an example let X (the extensive property) be Mass (M). 

In the transformation equation relating the Material Derivative 
of an extensive property (X)to the Control Volume, in 
Equation *1.28, X is mass, hence x is mass per unit mass, 
or X = 1. Thus Equation *1,28 becomes: 



dM 
dt _ 


9 

at 


> 

Q. 


C • S • 


pV.dA 


dt 


= 0 

material 

derivative 








SO , 

0 


a 

at 


/// . .. 
c,v. 


c . s . 


pV.dA 



This is the CONTINUITY EQUATION for a control volume, 
are other forms of this equation for special conditions 
as steady flov/, etc. which we will now develop. 



There 

such 






( 

i 

I 
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1.3 Control Mass versus Control Volume 



COMMENT - The concept developed In 1.2 is presented in 
an abbreviated form by the Text. Compare these tvjo 
parallel developments as this important concept must 
be thoroughly understood. 

a) Read section 1.6 of the Basic text. 

b) Can you interpret the conversion from the consideration 
of a fixed MASS system and the properties associated 
with it to a CONTROL VOLUME approach? VJhich is more 
convenient and WHY? 

c) Do you understand the significance of each term in 
Equation (1.7)?^ What does each term mean (Explain in 
words ) . 

d) Shov; v/hich term will drop out of Equation (1.7) in 
any steady flow situation and explain why. 



Note 1: Equation numbers v/ithout an asterisk refer to ' 

equations in the text (JOHN). 

Equations with an asterisk (e.g., *1.23) refer 

to equations in the study guides. 
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lA CONSERVATION OF MASS 



a. ) Read Section 1.7 in the text, 

COMMENT - In this course we shall assume that Matter 
and Energy are each separately conserved, as we 
shall not consider transformations such as occur 
in nuclear reactions. 

b. ) What is the Principal of Conservation of Mass, in 

words ? 

c. ) Now skip ahead and read Section 1.12. The real 

importance of the One-Dimensional Plov/ assumption 
is that it permits us to easily obtain approximate 
solutions (with acceptable engineering accuracy) 
to problems whose exact solutions would be extremely 
difficult. Any problem solution based on a simpli- 
fying assumption will yield results which vary to 
some extent from the real-world problem. The decision 
to utilize a simplification depends on its correlation 
to experimental results, the availability and 
cost of alternate methods of solution and the 
accuracy required of the solution. 

d. ) Now return to pg. 8 of the text and determine in 

detail the restrictions imposed in progressing from 
Eq. (1. 9) to Eq. (1.10) , 
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NOTE VJELL: The condition of ONE-DIMENSIONAL Flow 



imposes ^ restrictions on density, area, velocity, 
(or any other fluid property) from section to 
section in the direction of flov;, but does say that 
for a given location at a given instant of time, 
the flov; variables are assumed constant over 
the entire cross-section. 

e. ) Equation (1.10) is often written pAV = constant, 

in which case the V is assumed to be the component 
of velocity perpendicular to the area. V/hat quantity 
is the constant? 

HINT; Work out the Units 

f . ) Is it correct to say that the continuity equation 

can be stated as "The flow rate into a control 
volume equals the flow rate out of the control 
volume?" Is this mass flov; rate or volumetric flow 
rate? What condition must exist before the above 
statement is true? Pollov; through examples 1.1, 

1.2. Can you LIST the assumptions made for each 
of these problems? 

g. ) It is more convenient for the case of multiple 
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to express the continuity equation for one-dimen- 
sional, steady flow as: ^pAV = 0 
h.) An alternate form of the continuity equation (for 
one-dimensional steady flow) is obtained by first 
taking the natural logarithm of [pAV = constant] 
and then differentiating the result. Try it I 
You should get : 

^ ^ ^ = 0 (*l.ill) 

p V A 

This form is useful in interpreting the changes 
that must occur as fluid flows through a duct, 
channel or stream tube. It indicates that if 
mass is to be conserved, the changes in density, 
velocity and cross-sectional area must compensate 
for one another. For example, if the area is 
constant (dA = 0), then any increase in velocity 
must be accompanied by a corresponding decrease 
in density. This form of the continuity equation 
will also be useful in future derivations. 



1.5 Conservation of Energy 



a. ) Read Section 1.9 of the text. 

Section 1.8 v;lll be considered after the ENERGY con- 
cepts are firm. 

b. ) Do you recall the essential difference between a 

PATH FUNCTION and a STATE FUNCTION? The key lies 
In the definition of PATH functions which are only 
defined In crossing a system boundary . "Heat " and 
"pressure" are examples respectively of a PATH 
and a STATE function. List as many of each as you 
can recall. 

PATH FUNCTIONS STATE FUNCTIONS 



c.) What notation Is used to represent the Infini- 
tesimal change In a PATH function by JOHN? How 
about the notation for an Infinitesimal change In 
a STATE FUNCTION? V/hy the difference? 

Example: 6Q = 6W + dE 

Note that In equations l.l6 and subsequent the 
expression dQ/dt uses the same differential symbol 
for Heat (Q), a PATH function and for Energy (E) 
a STATE function. This is valid because we are 
now dealing v/ith a control volume, so the change 
In Heat (or work) is now a function of time only. 
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d. ) Two important concepts included in the development 

of equation (1,19) are; 

1 ) PLOW WORK 

2 ) ENTHALPY 

Flow work is the work done to push the fluid into 
(or out of) the system (i.e., CONTROL VOLUME). 

Since it is often difficult to measure Plov/ work 
being done on the system (v/hereas the Shaft work 
and other types of v.’ork can usually be measured), 
it is advantageous to combine Plov; work with specific 
Internal energy (a thermodynamic property) and 
thus define ENTHALPY (h). 

e. ) The formulation of the general ENERGY Equation as 

expressed in equation (1.19) can often be simpli- 
fied to facilitate obtaining a solution to practical 
problems in engineering. Consider the case of 
STEADY PLOV/. We can summarize the effects of 
STEADY PLOW by two conditions: 

1) MASS PLOW IN equals the MASS FLOW OUT of 
the CONTROL VOLUME. Why is this so? Does 
this define steady flov/? 

2) a(anything) = 0 

8t 

This assumption causes the first term of eq. 1.19 
to vanish. Similarly, as in examples 1.5 and 1.6, 
real problems can be readily treated by judicious 
definition of the control volume, v;hich makes the 
surface integral in eq . (1.19) extremely easy to 



Note carefully the use 
of partial versus the 
ordinary derivative 



evaluate , 
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f . ) The general ENERGY equation may be viewed as 

analogous to a bank account. Each of the terms 
represents a certain type of transaction, v/hlch 
must be balanced against each other to balance 
the books In accordance with the Law of Conserva- 
tion of ENERGY. To consolidate your understanding 
of the effect of each term, next to each of the 
statements below, write the term from the GENERAL 
ENERGY EQUATION, Equation (1.19), whose function 
Is described. (-^ You may have to separate certain 
terms ) 

1) The time rate of change of total energy 

within the control volume 

2) The net efflux of 

KINETIC ENERGY 

3) The net efflux of POTENTIAL ENERGY 



associated with the earth's 
gravitational field. 



4 ) 


The 


net efflux of flov; 


work 




and 


Internal energy 




5 ) 


The 


rate at which heat 


enters 




the 


control volume 





6) The rate at which work (other than flow work) 
leaves the control volume 
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g) From Thermodynamic considerations entropy can be broken 
down as follows; 

ds H ds^ + dSj^ ('^l.Sl) 

where;) "ds^" represents entropy change caused by 
( External Heat Transfer (or HEAT crossing 
( the boundary). It is assumed that this Heat 
I is transferred Reversibly . 

) "ds^" represents the Entropy change caused 
I by all the Internal IRREVERSIBLE effects. 

THUS 1) dSg = -Tjr (This can be + or - depending 

on the sign of SQ) 

2) ds^ is always positive! 

Now recall the "INEQUALITY of CLAUSIUS." Develop this 
from 5^1.51. 




From your Reviev; Questions you can recall 
dSg = 0 Implies ADIABATIC 
ds^ = 0 implies REVERSIBLE 

and ds = 0 means ISENTROPIC , therefore if a situation 
both ADIABATIC and Reversible it is ISENTROPIC. 

h) To simplify the application of the energy equation (1.19) 
to many engineering applications we v;ill again apply the 
assumption of ONE-DIMENSIONAL FLOW. 
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Now for ONE-DIMENSIONAL steady flov; a differential 
control volume may be conceptualized as the dotted area 
in Figure *1.51. 




Figure *1.51 



Thus we may apply the Energy equation (1.19) to this 
control volume, and divide out the flow rate with the 
following result: 



6q = [ (p+dp) (v+dv) - pv] + 

+ [(u+du) + + si|±d^] 






(n.53) 



If v/e expand this and neglect terms of higher than linear 
order we now have the Energy Equation in differential form: 



4q = d(pv) + + du + 



(n.5^) 



Fill in the intermediate steps for exercise. 
Hov; would you define 6q and 6v;? 
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1) The stage is now set to develop the "PRESSURE ENERGY 



EQUATION . " From the Thermodynamic relation 

T ds = dh - vdp (*1.55) 

we can nov; substitute Equation *1.51 to obtain 

dh = Tds^ + Tds. + ^ (*1.56) 

6 Ip 

Novf by substituting from equation *1.56 for dh in place 
of du + d(Pv) in equation *1.5^ and recognizing that 
Tds^ = 5q, we can obtain the useful form of the energy 
equation often called the PRESSURE ENERGY equation. 



+ yv_i + + Tds 



2g, 



Sr 



i ^^^shaft ° 



(*1.57) 



j) The PRESSURE ENERGY EQUATION is the source of a familiar 
relation from FLUID MECHANICS . If X'je consider the flow 
as REVERSIBLE (ds^^ = 0) and No shaft work crosses the 
boundary = 0) and if the fluid is considered 

incompressible (p =(|!! ) then we can integrate the Pressure 
Energy Equation between two points and obtain BERNOULLI ' s 
Equation. 



2g, 



gz- 



V 

2g, 



gz- 



(* 1 . 58 ) 



Summarize the assumptions made for Bernoulli's Equation 
for your own reviev;. 



1.6 - MOMENTUM 

a. Prior to the development of the equatlon(s) describing 
momentum, consider the concept of momentum — could you 
describe momentum to a 'non-engineering type'? Reflect 
on your own understanding and if you feel uncertain 
about this concept, break out your Physics text and 
refresh your understanding of 'MOMENTUM' and 'IMPULSE'. 

b. Study Section 1.8 in the text. Including examples 1.3 
and 1.4. Are we considering angular momentum? 

How about acceleration effects? 

c. Since neither angular momentum nor any acceleration 
effects are considered, we find that, from the text, 
equation (1.13) is in reality valid only under quite 
restricted conditions. Explain in your own v;ords the 
restriction related to each of the below listed concepts. 

a) Inertial Reference 

b) No Angular Momentum 

c) Control Volume (Fixed or Constant Velocity) 

Can you take equation (1.13) and "cross" both sides 
with the radius and come up with a similar relationship 
to deal with angular momentum? 

d. In Section 1.8 of the text, the material derivative of 
the extensive property momentum is being transformed to 
a control volume consideration, which is then related 
to the forces acting on the control volume by Newton's 
Second Law of Motion. VJe can apply to equation (1.13) 
in the text, the simplifying assumption of Steady and 
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One-dimensional flow. For One-dlmenslonal flow, the second 
term on the right side 




which describes the net amount of momentum per unit mass 
crossing the control boundary simplifies as follows: 



// 

: . s . 



p(V.dA) 

^c 





in 



(* 1 . 6 , 1 ) 



Prom Section 1.5 we can recall the two useful effects of 
Steady FIovj: 



1^ (anything) = 0 



and : 






Mass flow in 




flow out 




m 
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Therefore Equation (1.13) becomes 




fluid 






(n.6,2) 



e. Having obtained a relation which is now expressed in 

terms of the mass flow rate and velocities in and out of 
the control volume, we need now to examine the Forces 
acting on the control volume. 



FLOW 











w 



Figure *1.6-2 



In Figure *1.6-2 v;e see summarized the various forces 
acting on the control volume, where: 

W = Weight of Fluid in the control volume 
Fp = Upstream pressure forces 

F = Downstream pressure forces 

P2 

F = Pressure forces exerted by the walls 
Pw 

F„ = Friction forces exerted by the ivall 
^ w 

It is often common practice to group the last two together 
and refer to them as 'Enclosure Forces'. VJhere necessary 
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we can also include other forces such as electrical, 
magnetic or other types of field forces which for 
simplicity we will not consider here. 

V/e establish co-ordinates and designate the various 
forces by subscript acting on a differential control 
volume of length dx. 




Figure ^1.6-3 

Since this is steady, ONE-DIMENSIONAL flow in the 
x-direction, the sum of Forces in any other direction 
must equal zero from equation (*1.6-2). 
f. We can novj summarize the x-components of the Forces. 

Note we vrill establish the sign convention that force 
components in the upstream direction are positive . 
Upstream Pressure Forces 

F, = PA (*1.6,3) 

lx 



5 ^ 



Downstream Pressure Forces 



^2x = - (P + dP)(A + dA) (*1.6,10 

Expanding and Neglecting higher ordered terms: 

Ppx = - [PA + PdA + AdP] -(*1.6,5) 

Wall Pressure Forces 

F 3 X = (p + ^)(^dL) sin e (* 1 . 6 , 6 ) 

where is the Mean Wetted Perimeter 
and 0 is the divergence angle of the wall 
if we recognize that 

( l^dL) sin 0 = dA (*1.6,7) 

and we expand and neglect higher ordered terms. 

F, = PdA (*1.6,8) 

Friction Forces (J) 

If we define as the mean shear stress along 

the wall 

Fi^x = - (fdL) cos 0 ] (*1.6,9) 

recognizing that dL cos 0 = dx 

Fi, = - T, ^ dx (*1.6,10) 

^x f ' 
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Gravity Force 



= - [ (A + dx ] [p + ^ cos (j) (*1.6,11) 

where <(> is the angle of inclination of the x axis 
but dx cos <j) = dz. Expanding and neglecting higher 
order terms: 



5x 



- Ap — dz 
^c 



(* 1 . 6 , 12 ) 



g. Having shown that the right side of (*1.6,2) is: 



— [V , - V. ] 

gp out in 



(*1.6, 2a) 



The X component of this becomes 



m 

S. 



[V + dV - V] 



(*1.6,2b) 



or 



m 



dV 



You should equate the forces to the right side, cancel 
terms and obtain the following result: 






AVdV 



AdP + T. F dx + p A dz + p = 0 

X 'CP 



g. 



(*1.6,13) 



To evaluate the Shear term we will introduce a non- 
dimensional friction co-efficient (f) which will relate 
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shear force to the dynamic pressure. (Sometimes the 
friction factor is defined v/ithout the factor of 4.) 



f E 



pv' 

2g. 



(* 1 . 6 , 1 ^) 



We also introduce the concept of an 'Equivalent diameter' 



E ^ 



Then 






(*1.6,15) 






(*1.6, 15a) 



and substituting these expressions into Equation (*1.6,13), 
divide the equation by density and we obtain a useful 
form of the Momentum Equation in the direction of fluid 
flow for steady, one-dimensional flow. 



dP 

p 



fV" 



Dp2g 



dx 



E &c 



+ dz 

Sc 



VdV 



= 0 



(* 1 . 6 , 16 ) 
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1.7 - THERMODYNAMIC CONSIDERATIONS and EQUATION of STATE 



a. Read Section 1.10. From your basic Thermodynamics course 
you may recall that entropy changes can fall into one of 
two different categories: 

ds = ds. + ds^ 

1 e 

where ds^^ is the change in entropy due to irreversible 
effects (such as friction) 

and dSg ( ^ ) is the change in entropy due to heat transfer. 
Comment on ds^^ 

dSj_ is ali-^ays POSITIVE since all Irreversible effects 
generate entropy. But if, in a process, the changes 
in pressure temperature, etc. are vanishingly small, 
then this process may be treated as reversible, 
hence dSj^ = 0. 

Comment on ds 

e 

dSg can be positive or negative (depending on the 
direction of the heat transfer) and is zero if no 
heat transfer occurs. 

b. Read Section 1.11. This section reviews some vital 
concepts and relations which will be vital in many 
applications. Example 1.7 provides an excellent opportu- 
nity to follow through a problem using a systematic 
problem solution technique. 
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A recommended method for approaching a problem is provided 
below. Although each of these factors may not figure in 
every problem, as was the case in example 1.7, they should 
help you organize a problem solution. 

1) Identify (and sketch) CONTROL VOLUME, the dotted 
area in Example 1.7. 

2) Identify FORCES acting on the matter inside the 
CONTROL VOLUME. 

3) Identify where the fluid enters and leaves the 
CONTROL VOLUME. 

4) Identify any energy (Q and W) being transferred 
. across the control surface. 

5) ESTABLISH A CO-ORDINATE SYSTEM. DON'T FORGET SIGN 
CONVENTIONS. 

6) BE ESPECIALLY CAREFUL WITH THE SIGNS OF VECTOR 
QUANTITIES SUCH AS V and F. 

7) EVALUATE WHAT IS KNOWN and WHAT IS UNKNOWN. DECIDE 
WHAT RELATIONSHIPS WILL BE USEFUL. 



The 


Basic Relationships v/ith 


which v;e will analyze 


Gas 


Dynamics will be 


: 






1) 


CONTINUITY 


(Derived 


from Conservation of Mass) 




2) 


ENERGY 


(Derived 


from Conservation of Energy) 




3) 


MOMENTUM 


(Derived 


from NEWTON'S SECOND LAW) 






An EQUATION 


of STATE 


p = (^ 










p = pRT 










others 




5) 


THERMODYNAMIC CONSIDERATIONS 



c. Revlev; Section 1.13. 
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UNIT 2 - STAGNATION CONCEPTS WAVE PROPAGATION and MACH NUMBERS 



- OBJECTIVES - 



After successfully completing this unit you should be able to: 

1. Explain the stagnation state concept, hov/ it is achieved 
and utilized. 

2. Compare the concept of a STAGNATION property vjith that of 
a STATIC property. 

3. Dravj a T - S diagram representing a flow system and indicate 
static and stagnation points for an arbitrary section. 

4. Define STAGNATION ENTHALPY. 

5. Introduce the stagnation concept into the "Pressure Energy" 
equation and derive the "Stagnation Pressure Energy" 
equation . 

6. Explain how sound is propagated through any medium (solid, 
liquid or gas). 

7. Define SONIC VELOCITY. 

8. Utilize a control volume analysis to derive the general 
expression for the velocity of wave propagation in an 
arbitrary medium, starting with the "Continuity" and 
"Momentum" equations for steady, one dimensional flov/. 

9. Recall the relations for: 

a. Speed of Sound in an arbitrary medium. 

b. Speed of Sound in a Perfect Gas. 

c. Mach Number. 

10. Discuss the propagation of signal waves from a moving body 
in a fluid, and explain what is meant by 'Zone of Action', 
'Zone of Silence', 'Mach Cone', and 'Mach Angle'. 

11. Graphically portray the simplified flow over a wedge 
shaped body for SUBSONIC and SUPERSONIC flow. 

12. Express the basic equations in terms of Mach number for 
a Perfect Gas by shovjing that : 

a. the "Continuity" equation can be written as 
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b. the "Energy" equation can be v/rltten as 



(1 


+ M^)dh 


+ hdM^ 


= 


6 q - 6 V7 


the 


"Momentum" 


equation can 


be 


vjritten as 


dP 

P 


+ J dM^ + 


X — + 

2 T 


Y 

2 


f'd X 

M - 0 



6l 






!0 
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UNIT 2 - STUDY GUIDE 



2 . 1 The Stagnation Concept 

a. V/e introduce the concept of a reference state defined as 
that thermodynamic state which would exist if the fluid 
v/ere brought to zero velocity and zero potential. To 
yield a consistent reference state we must qualify how 
this "stagnation process" would be accomplished. The 
stagnation state m^ust be reached: 

(a) without any energy exchange (Q = V/ = 0) 

• (b) without losses. 

By virtue of (a), dSe = 0, and from (b), dS^ = 0. Thus 
the stagnation process is isentropic. 

b. Vie can imagine the following example of actually carrying 
out the stagnation process. Consider fluid v/hich is 
flowing under condition 1 in Figure *2.1 (these condi- 
tions are referred to as the "static" conditions). At 

the fluid has been brought to zero velocity and zero 
potential under the above restrictions. 




c. Apply the "energy" equation to the streamtube and you have: 

0 0 0 0 

" si; " -iT "X-2 “ "2 "ifc " A " a-2 



which simplifies to: 



Vt gz. 



h, + ■ 5 — 
1 2 g 



c Sc 



= h. 



(* 2 . 1 , 2 ) 



But condition represents the "sta^ation state" 

corresponding to the "static state" • We thus call 
hp the s tagnation or total enthalpy c^respondlng to 
state and designate it as h^ . 
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Thus : 



h 



t 



1 





gz- 



g. 



(* 2 . 1 , 3 ) 



Or, for any static state we have in general: 

.2 



= h + 



V 

2g. 






(* 2 . 1 , 5 ) 



The introduction of the stagnation (or total) enthalpy 
makes it possible to write equations in a simplified 
form. For example, the steady flow energy equation 
becomes : 



dq = dWg + dh^ 



(* 2 . 1 , 6 ) 



or 



+ q 



1-2 




+ 




(* 2 . 1 , 7 ) 



You should note that the stagnation state is a reference 
state and may or may not actually exist in the flow 
system. Also, in general, each point in a flow system 
has a different stagnation state as shown below: 



T 




© 



Physical System 




S 



Figure ^2.2 



! 

I 

\ 



( 



Also one must realize that when the frame of reference 
is changed then stagnation conditions change, although 
the static conditions remain the same. (Static properties 
are defined as those that would be measured if the 
measuring devices move v;ith the fluid.) 

Consider still air and the earth as a reference frame. 



V = 0 P = 1^4.7 psia T = 520°R 




Figure *2.3 



In this case, since the velocity is zero (with respect to 
the frame of reference), the static and stagnation 
conditions are the same. 



f. Nov; let’s change the frame of reference by flying through 
the air on a missile at 600 ft/sec. As v;e look forward it 
appears that the air is coming at us at 600 ft/sec. Its 
static pressure and temperature are unchanged at 1^.7 psia 
and 520° R respectively. 



V = 600 ft/sec 
P = 1^4.7 psia 5 ^ P^ 
T = 520° B. ^ 




But now the air has a velocity (v;ith respect to the frame 
of reference) and thus the static and stagnation condi- 
tions are different. 

g. You will soon learn how to compute the stagnation condi- 
tions. Incidentally, is there any place in this last 
system where the stagnation conditions actually exist? 
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2 . 2 The Stagnation Pressure Energy Equation 



a. Consider two stations in a flov; system that are very 

close together and with thermodynamic states differentially 
separated as shovjn in Figure *2.5. Also shown are the 
corresponding stagnation states for these two stations. 




Figure *2.5 



We may write the following property relation between 
points 1 and 2: 



Tds = dh - vdp 



(* 2 . 2 , 1 ) 



Since this is a valid property relation we may also write 
this between points 1^ and 2^. 



ct 

a> 

Ui 

ct 

II 


dh^ - 




(*2.2,2) 


However , 




ds 




and 


ds = 


dSe + 




Thus we may 


v;rite 


• 






dSj^) 


' - ''t'^Pt 


(*2.2,3) 



Recall the "energy" equation written in the form; 

8q = 8w + dh, (*2.2,4) 

S 1/ 
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Substitute equation *2.2,3 into equation *2.2,4 and you 
should obtain: 



5q = 6v;g + T^(dSg + ds^) + v^dp^ (*2.2,5) 

Also recall that: 6q = Tds^ (*2.2,6) 

By substituting equation *2.2,6 into *2.2,5 and noting 
that — 




you should obtain the following equation which is called 
the "Stagnation Pressure Energy" equation: 



— + 5w 
t 



+ ds (T.-T) + T.ds. = 0 

e t t 1 



(«2.2,7) 



b. Consider what happens v;hen 



(a) 


There 


is 


no 


work transfer 


8Wg 


= 0 


(b) 


There 


is 


no 


heat transfer 


dse 


= 0 


(c) 


There 


are 


nc 


1 losses 


ds^ 


= 0 



Under these conditions = 0 

pt 



or dP^ = 0 (*2.2,8) 

or P^ = constant (*2.2,9) 



Note that in general the total pressure will not remain 
constant; only under a special set of circumstances will 
equation *2.2,9 hold true. 



66 



2 . 3 Velocity of Sound 



a. Read Sections 2.1, 2.2 and 2.3 in the text. 

b. Recall from your previous v/ork in fluid mechanics, the 
definition of a stream line. V/hat, if any, properties 
reamin constant along a stream line? V/hat significance 
does the distance betv/een stream lines have in an 
Incompressible medium? 

c. Consider carefully the implications of section 2.2 with 
regard to the velocity of sound in a compressible medium. 

It is important to note that a sound wave results from an 
infinitesimally small impulse. The text uses a moving 
control volume v.'hich encloses the wave to analyze this 
unsteady situation as a steady flow problem. We can view 
this process as a simple change of reference from that of 
a fixed observer viewing the moving wave to that of an 
observer v;ho is "riding" the wave. The moving observer 
experiences the fluid flowing from right to left in 
Figure 2.5 at a velocity of sound (a) upstream and 
velocity (V-a) downstream. This does not change the 
static conditions. V/hat about the stagnation conditions? 

d. Follow through Example 2.1. The author of our text 
introduces the relationship of 'compressibilities' which 
may be unfamiliar to you. V/hat other ratio is expressed 
by (y)? We need not concern ourselves at this stage with 
situations in which air does not behave approximately as 
a perfect gas, other than to examine each problem with 
this in mind. In general, for flow situations vjhere 
extremely high pressures and/or extremely low temperatures 
are encountered the 'Perfect Gas' approximation is not 
valid. Otherwise v/e may treat air and other gasses (if 

no appreciable molecular interaction is taking place) as 
a perfect gas. 

e. Read sections 2.^4 and 2.5 in the text. In Figure 2.10 in 
the text, the circles indicate the location of the original 
sound wave emitted while at point 'O'. Figure *2.6 in the 
Study Guide shows the location of the sound v.^ave initially 
emitted at times, t = 0, 1, & 2 as they appear at time, 

t = 3. 
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f. Notice the distance betv;een successive sound v;aves is 
less in the direction of travel of the point projectile. 
Now compare Figure ^ 2.6 to Figure 2.11 in the text. 

In Figure 2.11 the point is overtaking the sound waves 
it emits as it moves. The v^^aves are left behind creating 
a situation in a sense like the wave created by dragging 
the point of a stick rapidly across the surface of a 
still body of v;ater. The Mach Wave (or Mach Cone) 
produced by a point source is an infinitely vjeak pressure 
wave, but the addition of a multitude of such waves 
produced by a finite body results in the shock shown in 
Figure 2.13 in the text. 

g. It should be emphasized again that both velocity and 
speed of sound are local conditions. 



Now consider the problem below: 



air 




Note: No information 

may be assumed about- 
pressure, area, 
temperature, etc.- 
changes from 1 to 2. 



= 500 ft/sec = 1000 ft/sec 



M^ = .5 



M = ? 
2 



Therefore : 



a^ = 1000 ft/sec 



a2 = ? 



Is enough information given to compute a 2 ? 



68 



t ■...• 






11 ;^ 










■ 







J 





{ 



2.k - Basic Equations In Terms of Mach Number 



a. The Importance of Mach Number as a parameter vrhlch 

controls the behavior of fluid flovj suggests the value 
of expressing the basic equations developed In Unit 1 
In terms of Mach Number as v/ell as velocity. This can 



be done 
for the 


for steady, one dimensional flow quite readily, 
case of perfect gases. 


b. For the 


CONTINUITY equation — 



Equation (1.10) from the text may be expressed for steady 
flow as : 

m = pAV = (f (*2.il,l) 

where (m) Is mass flow rate. 

From the Perfect Gas Lav;: 



P 


= ^ («2.4,2) 


and the 


definition of Mach Number: 


V 


= Ma (*2.4,3) 



Now recall the expression for sonic velocity In a Perfect 



Gas : 




a 


= (*2.4,4) 


where 


Y, g and R are defined In Section 1.11 of 
the text. 


We thus 


have : 


• 

m = 


pAV = = PAM*^-^ (*2.4,5) 



Thus, for the steady, one dimensional flow of a Perfect 
Gas, the "Continuity" equation becomes: 

= ({: (*2.il,6) 
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c. Differentiating the "Energy" equation. Equation (1.19) in 
the text, and neglecting the potential term, v;e can show 
for steady flov; that; 



d[h + 


2g 3 - 6q - <5Wg (*2.4,8) 


You saw 


that for a perfect gas: 


a2 = 


Yg^RT (^2.4,9) 



hence 



d[h + 


^2 “ 5q - <5Wg (*2.4,10) 



We can re-write this expression as: 

d[h + ( ^ ) = 6q - 6Wg (*2.i|,ll) 

From section 1.11 recall that for a perfect gas: 



1 — 1 

1 


CM 
1 — 1 

CM 

II 



and 



dh = 


c dT (*2.4,13) 



Show that 



d[h(l 


+ ^ M^)] = 6q - 6Wg (*2.4,14) 


Now use 


the product rule for differentials to show that: 



(1 + ^ M^) dh + ^ hdM^ = 6q - Sw^ (*2.4,15) 
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d. We derived in Unit 1.6 of the Study Guide an expression 
for the "Momentum” equation in steady, one dimensional 
flow : 



P 



+ 



fV^ 



dx 



Sc 



0 



(*2. 4,16) 



Divide through by P, multiply each term by p, neglect the 
dz term, and rev;rlte the last term to yield: 



dP pfV^ pdV^ 

P PDgPgc P2g^ 



(*2.4,17) 



Now, substitute pRT for P, multiply the second and third 
terms by y/y and regroup the second term to obtain: 



P 



+ 





+ YdV^ 
2yg^RT 



0 



(*2.4,18) 



The third term can be modified by considering the product 
rules for the differential of the tvjo variables M and T: 



Y dV^ 



2 



Y d(a^M^) 
2 a2 




+ J ^ 



(*2.4,19) 



So now the "Momentum" equation may be v;ritten: 



P 



+ 




X 

2 



M 



2 ^ 
T 



I M^(f 1^) 



0 (*2.4,20) 
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UNIT 3 - ISENTROPIC FLOV/ - OBJECTIVES 



The student shall be able to: 

1. Define an ISENTROPIC PROCESS and explain the relation- 
ship among reversible, adiabatic and isentropic processes. 

2. Show graphically how pressure, density and area vary in 
steady, one-dimensional, isentropic flow as Mach Number 
ranges from zero to supersonic values. 

3. Explain the difference between STATIC and STAGNATION 
properties. Given all the static properties, compute 
the stagnation properties. 

Describe v;hat is meant by a "choked" flov; passage. 

5. Compare the function of a NOZZLE and a DIFFUSER. Sketch 
physical devices that perform as each for subsonic 

and supersonic flow. 

6. Simplify the basic equations for CONTINUITY, ENERGY 
and MOMENTUM to relate differential changes in density, 
pressure and volum.e to a differential change in area 
for steady, one-dimensional flow through a varying area 
passage . 

7. Define the reference condition, (*) and the properties 
associated v;ith it. (l.e.. A*, P*, T* , p*, etc.) 

8. Define stagnation temperature (T^) of a perfect gas in 
terms of temperature (T) , Mach Number (M) and ratio of 
specific heats (y)* 

9. Define stagnation pressure (Pt ) of a perfect gas in 
terms of temperature (T), Mach Number (M) and ratio of 
specific heats (y). 

10. Define stagnation enthalpy (ht ) in terms of enthalpy (h) 
and velocity (V). Put in term.s of M and y for a perfect 
gas . 

11. Depict graphically the general relationships of velocity, 
density and area in varying area adiabatic flow where 
the Mach Number varies from zero to supersonic values. 

12. Express the loss relationship in adiabatic flow (ASj_) as 
a function of stagnation pressures (Pt) or reference 
areas (A*) between two points in the flov;. 
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13 . Cite the necessary conditions for the reference area 

A* to remain constant for all points in the flov; field. 

1 4 . Derive the working equations for a perfect gas relating 
property ratios betv^een two points, in both adiabatic 
and isentropic flov/s, as a function of Mach Number (M) , 
ratio of specific heats (y) and change in entropy (As). 

15. State and interpret the relation between stagnation 
pressure (Pt) and the reference area (A*) for the 
process between two points in adiabatic flov/. 

16. Utilize the adiabatic and isentropic flow relations 

and the isentropic tables to solve typical flow problems. 
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STUDY GUIDE - UNIT 3 



ADIABATIC AND ISENTROPIC FLOW 



\ \ 


\ 


Engineering problems typically Involve many factors which 
affect the fluid flov/ conditions and frequently one of 
these factors is predominant. In this case it is 
possible to consider only the controlling factor and 
develop a simple solution to the problem. These 
approximate solutions can be applied with great confidence 
to many engineering situations. To the engineering 
student, this approach also provides insight concerning 
the significance, magnitude and relative Importance 
of each of the factors considered, as the basic state 
variables are changed. 



3.1 - Variable Area Flow In General 



a. We will now consider flow in a variable area channel 
without heat transfer, or "adiabatic" flov;. Many real 
life situations fit into this category. We have 
learned that for many situations air and other gases 
flowing in a continuous field obey the Perfect Gas Lai*; 
reasonably well. Since this is the case, we will 
develop v;orking equations using the 'perfect gas' 
assumption. Additionally, we are motivated, for the 
sake of simplicity, to continue to assume steady, one- 
dimensional flow with no work (that is, shaft work) 
done on or by the system. 

If we consider the adiabatic flow as it occurs in real 
situations, some accounting must be made for losses 
which occur as a result of friction and other irreversi- 
ble factors. When these factors are considered negli- 
gible and the flov; is taken to be reversible, the flov; 
meets the criteria for isentropic flov/. This concept 
was covered in brief in the Review Questions. We may 
viev; the Isentropic flov; case as the 'ideal' or 
'Standard' with which to measure all adiabatic flov/. 

b. Read sections 3*1 and 3*2 in the text. 

c. What common simplification is used to obtain equations 
3.1 and 3-2 from the basic equations? When is it 
justifiable to neglect the higher ordered terms and 
assume that these versions of the basic equations 
(3.1, 3.2, etc.) are correct? 
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d. To gain some insight into the manner in v/hich the flow 

variables are changing v;hen a fluid encounters a variable 
area, it v/ill be convenient to restrict the analysis to 
flow with no losses. V7e shall start with the energy 
equation : 



8.q = 8Wg + dh^ (*3.1,1) 

Since we have made the assumptions of no heat transfer 
or shaft v;ork : 



dhj = 0 



(* 3 . 1 , 2 ) 



Recalling the definition of stagnation enthalpy (neglecting 
potential energy change): 



h 



t 



h + 




Differentiating: 



(* 3 . 1 , 3 ) 



dh^ 



dh + 



VdV 



but : 



dhj = 0 



so : 



dh 



VdV 

Sc 



(* 3 . 1 ,^) 



Now, we may start with equation (1.22), from the Second 
Law of Thermodynamics: 



Tds 



dh 



6? 

P 



(* 3 . 1 , 5 ) 



Since this process is both reversible and adiabatic, it 
is isentropic . Thus ds = 0. 
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So finally: 



dh = ^ (*3.1,6) 

Considering the equations (*3.1,^) and (*3.1,6) we obtain 

VdV ^ 

Sc P 

or 

(*3.1,7) 

equation 3.1 and the differential 
equation becomes; 



dV = - 



g^dP 

pV 



We introduce this into 
form of the continuity 



dp , dA g dP 

v;hich can be solved for 



0 



^ r Ip + M 1 

P Sc P A ^ 

The definition of sonic velocity; 



S. 



can be written 



3p 

for 



s=constant 
our case as : 



"■c dp 

because this flov; is isentropic 



Now, rearranging equation 2.3: 

.2 

dP = — dp 

Sc 



(* 3 . 1 , 8 ) 



(*3.1,9) 



(2.3) 



(* 3 . 1 , 10 ) 
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Substituting this expression for dP into equation *3.1,9 
yields : 



d£ 

p 



Y r I dA 1 

a2 ^ 0 * 



(* 3 . 1 , 11 ) 



From the definition of Mach Number 




show that : 



ip 

p 



[ 




] 



A 



(«3.1>12) 



If we now substitute equation *3.1,12 into equation 3-1 
in the the text you can shov/ that: 



6N 

V 



1-M 



1 M 

2 ^ A 



But we recall that: 
~®c 

dV = — dP 

pV 



(*3.1>13) 



(*3.1,7) 



or 

^ dP (*3.1,7a) 

pV 

d. A/" 

Substitute this expression for into equation *3.1,13 
and solve for dP: 



dP 



pV^ r 1 T ^ 
®c 1-M^ ^ 



(3.5) 
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We now collect the principal relations which will be 
referred to in the follovjing sections. 



dP 



pvf r 1 1 ^ 

®c 1-M^ ^ 



(3.5) 



dp 

p 




1-M 



2 




(* 3 . 1 , 12 ) 



V 



1 

2 

1-M"^ 




(*3.1,13) 



e. Consider now what must be happening when fluid flows 
through the variable area channel. Fill in the blanks 
below. Assume that the pressure is always decreasing. 
Thus dP is negative. 

From equation 3*5 you see that if M< 1, dA must be 

indicating that the area is , whereas if M> 1, 

dA must be and the area is . 

From here you can move to equation *3.1,12. If M< 1 and 



dA 


then 


dp 


must 


be 


. If M > 1 and 


dA 


then 


dp 


must 


be 





Looking at equation *3.1,13 reveals that if M< 1 and 

dA then dV must be meaning that velocity 

is , vjhereas if M > 1 and dA then dV 

must be and velocity is . 
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f. We can summarize the effects and variation of velocity, 
density and area as Mach Number varies from zero to 
supersonic values. Refer to Figure *3-1 and recall the 
continuity equation: 



^ ^ ^ = 0 (n.ni) 

p A V 




Summary of variation of p , V and A. 



M << 1 


M = 1 


M > 1 


1. Density^ constant 


1. Density 


1. Density 




compensates 


compensates 


2. Area 


Velocity 


Area 


compensates 




and 


Velocity 


2. dA = 0 


Velocity 



g. Now read sections 3-3 and 3-^ in the text. Consider 
carefully the results pictured in Figures 3-3j 3. ■4 and 
3 . 5 . The situation depicted in Figure 3-5 in the text 
can only exist if the proper pressure conditions exist 
throughout the nozzle. 
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h. A 'nozzle' is a device which converts enthalpy (or pressure 
energy, for the case of incompressible fluids) into 
kinetic energy. A 'diffuser' is a device vjhlch converts 
kinetic energy into enthalpy. This distinction is often 
overlooked and a device is labeled a "nozzle" simply 
because it is a convergent section, with no regard for 
the flov; conditions. 



3.2 - Computation of Stagnation Conditions for a Perfect Gas 



a. Read the first part of Section 3-4 (from page 39 to the 
middle of page 4l). 

b. Particularly note equations 3*6 and 3- 7* These will be 
used frequently and should be learned. 



3.3 - Adiabatic Flow for a Perfect Gas 



a. We now return to the problem of adiabatic varying area 
flow and include any losses that might exist. To make 
this problem tractable vie will assume that the fluid 
is a perfect gas. 

b. Prom earlier work recall the follov/ing relations which 
hold true for a perfect gas : 



Ah 



c AT 
P 



(1.25) 



c 



c 



R 



( 1 . 26 ) 



P 



V 



c 



p 



Y 



c 



V 



Prom section 
1.11 in text 



c 



P 




Prom 3.1 we showed that 



8q 



6Ws + dh^ 



(*3.1,1) 
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So for adiabatic flow with no v;ork : 



dh 



t 



0 



(* 3 . 1 , 2 ) 



or 



h 



t 



$ 



This is true for any fluid under conditions where Q = V/^ = 0. 
From equation 1.25, if one considers a perfect gas, 
the stagnation temperature (T^) is also constant. 



c. Review the relations we already know to determine which 
ones may be useful in analyzing adiabatic flow. First 
consider the continuity equation for steady, one- 
dimensional flow: 



m 



pAV 




(*3.3,1) 



Examine the flow between two arbitrary points in a flow 
system such as in Figure ^3.2. 



FLOW 



Figure *3.2 



1 



2 



Equation *3.3,1 may be written: 





(*3.3,2) 



or 



A 



A 



1 



2 




(*3.3,3) 
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From the definition of Mach Number: 



V = Ma 






and the expression for the speed of sound in a perfect 
gas : 



a = ^/yg^RT 



(*3.3,5) 



and recalling the equation of state for a perfect gas: 



P = pRT 



(*3.3,6) 



we can re-write equation *3.3,3 as: 




PlMi 

P2M2 




(*3.3,7) 



Work out the intermediate steps from *3.3,3 to *3.3,7 in 
order to corroborate this expression. 



d. We will now turn to the energy equation. We have seen 
that : 



h 



t 



h + 




(* 3 . 3 , 8 ) 



Hence the energy equation may be written between two 
points as: 



h 



t 



1 



+ 



^ 1-^2 




+ 




(* 3 . 3 , 9 ) 



Since in a perfect gas under adiabatic flovj 

conditions with no shaft work, you can now re-write equation 
1.25 from the text, betv;een any arbitrary point in the 
system and stagnation conditions. You should get: 
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i 

1 



h 



Cp [Tp - T] 



(*3.3,10) 



h*. - 



Now using equation *3*3j8 we can shov; that 



= T [1 + ^ M^] 



(3.6) 



The development of this equation is found in section 
3.^ in the text. 

Since is a constant betv/een any two points in the flow, 
v;e can now write: 



T. 



[1 + ^ = 



T 2 [1 + ^ 



or 



T. 



v-1 2 

1 + VMi 



1 + ^ 



(*3.3,11) 



(* 3 . 3 , 12 ) 



e. From the stagnation energy equation for adiabatic flovj 
of a perfect gas with no v;ork : 




ds . 
1 

R 



0 



(*3.3,13) 



If v;e integrate both sides between two points (1 & 2), 
we get : 



In 





R 




(* 3 . 3 , 1 ^) 



Since dsg = 0, Asi = As so, taking the anti-log of 
both sides of equation *3.3,1^ yields: 





As 


1 : 

C\J 

.p 

IP 


" R 


1 

1 — 1 
p 
p 





(*3.3,15) 



This expression is of particular value. 
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f. stagnation pressure is defined as the pressure obtained 
if the flow is brought to rest 1 sen tropic ally . 

Between two points we can recall from section 1 (for an 
Isentropic process): 




Let condition 2 be the stagnation state and condition 1 
be the static state. Use equations 1.30 and 3-6 to show 
that : 



^ = (1 + ^ ^ ( 3 . 7 ) 

Thus between two points: 



P 



t2 



P 



tl 



Ppd + ^ M2^)^“^ 
Ppd + ^ 




Rearrange this equation 




1 + ^ 
1 + ^ 



to shov: that : 




(•»3.3d6) 



(*3.3d7) 



We can now utilize equation *3>3d7 to obtain a relation 
which expresses the change in area in terms of Mach 
Number and a loss term. 

From continuity show that: 




P^Mi 

PpMp 




(*3.3d8) 
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You can now substitute from equations *3-3jl7 and *3.3,12 
to obtain: 



A 



2 



A 



1 



Y 1/2 



r 1 + 1 


Y-1 / 


M, 


" 1 + M ^ 1 


2 2 


(p ^ 1 


1 


2 1 


1 + 1^ M ^ 


V / 


M2 




2 1- 






2 2 -I 



(*3.3,19) 



By combining terms, equation *3-3,19 can be v;ritten; 





1 + 



L 1 + 



Y-1 

2 

Y-1 

2 




Y + 1 

2 (y- 1 ) / As 




(*3.3,20) 



3.4 - Isentroplc Tables 



d. Read the remainder of section 3.4 in the text. Note the 
introduction of a new reference condition denoted by *. 

■At this location the Mach Number is unity and this 
condition can be reached by many processes. This is a 
hypothetical reference condition v;hich may or may not 
exist in the system (as v/as the case for the stagnation 
condition ) . 

b. One can now simplify many general equations that have been 
developed by assuming no losses. 

Consider an isentropic flov; case viev/ed in the T — S plane. 
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By utilizing the expressions for area, pressure and 
temperature ratios and letting one state point be a 
reference condition, if vie knov; v;hat the fluid is (i.e., 
the value of y) then we can construct tables v;hich relate 
these ratios to various values of Mach Number for the 
isentropic case. 

As an example, let 2 be an arbitrary point in the flow 
system and let 1 be the stagnation reference state. Then 

T 2 = T (hence M 2 = M = any value) 

(hence M^ = 0) 

and equation *3*3jl2 can be written as: 




1 + M^ 



f(M,y) 



(* 3 .^, 1 ) 



As another example, let 2 be an arbitrary point in the 
flow system and let 1 be the reference state *. 

Then: A 2 = A (hence M 2 = M = any value) 

A^ = A* (hence M^^ = 1) 



Remembering that v;e are considering an isentropic process, 
equation ^3-3,20 can be written as: 



A 

W 



1 

M 



1 + M^ 

y+T 

2 



Y +1 

■2TY-iy 



f(M,y) 



(*3.ii,2) 



c. We see that for isentropic processes these Important 
ratios are simply functions of Mach Number and gamma (y) 
which lead to tabulation of them in Appendix A. 

d. Evaluate equation *3*^jl for an arbitrary value of M> 1 
and an arbitrary value of M<1. Compare your answers 
with the values for T/Tt in the appendix. Notice that 
the tables are dimensionless and are alv;ays expressed in 
terms of a ratio of pressures, temperatures, etc. 
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3.5 - Relations Betv/een Pt, and A* 



a. V/e can now relate the stagnation pressure (P^) and the 
sonic reference area (A^). Take Ap as A]^* and take 
A 2 as A 2 * in equation *3.3,20. By definition = Mp = 1. 
This expression reduces to 





(*3.5,1) 



Vie sav; that : 





(*3.3,15) 



Thus 




or 



^t2 ^2* ~ ^tl^l* ~ constant 



(*3.5,3) 



Note that equations *3.5,2 and *3.5,3 apply to any 
adiabatic situation of flow of a perfect gas where 
«s = o- 

b. It should be noted (and we can verify from equation *3.3,20) 
that for two different points in the flow field (1 & 2), 

A-]_* is only equal to A 2 * if the flow is isentropic betv;een 
points 1 & 2. 



c. Now read sections 3.5, 3.6 and 3-1 in the text. 
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UNIT i] - NORMAL and OBLIQUE SHOCKS - Objectives 



The student shall be able to: 

1. List the assumptions used to analyze a standing normal 
shock . 

2. Sketch a "shock process" on a T-S diagram, noting the 
pertinent features of the phenomenom. 

3. Explain why an "expansion shock v;ave" cannot exist. 

Recognize the necessary and sufficient conditions to 
solve "normal shock" problems. 

5. Derive the expression which describes the properties 
on the dovmstream side of a standing normal shock in 
terms of properties on the upstream side. 

6. Explain how shock tables may be developed from the 
relations describing the properties on the downstream 
side of a standing normal shock in terms of the properties 
on the upstream side. 

7. Demonstrate the ability to solve typical normal shock 
problems by use of the tables and/or normal shock 
equations. 

8. Explain how an oblique shock can be described by the 
superposition of a normal shock and another flow field. 

9. Identify which properties remain constant and which 
change vjhen a uniform velocity is superimposed on the 
flow field. 

10. Show (by diagrams) how the "shock angle" and the 
"deflection angle" are defined. 

11. Describe the general results of an oblique shock analysis 
in terms of a diagram such as "shock angle" versus 
"Miach number" for various deflection angles. 

12. Distinguish between weak and strong shocks. Knov; what 
conditions cause each to form. 

13 . Describe the conditions v.'hich cause a detached shock 
to form. 

iH . Solve typical problems involving oblique shocks (such 
as air inlets, nozzle outlets, wedges, etc.) 
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15. Describe the treatment of a moving normal shock wave 
so as to apply the relations developed for the 
standing normal shock. 






I 

4 

a 



r;- 
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Unit ^ - NORMAL and OBLIQUE SHOCKS 



U.l - The Shock Mechanism 



a. To this point we have treated fluid flow problems which 
occured in a continuous medium; that is, v/herein no discon- 
tinuities existed within the control volume under consider- 
ation. V/e now turn to consideration of a discontinuous 
process - the "shock". A shock v.'ave characteristically 
occurs v/ithin a very thin, but finite volume. The thick- 
ness of a shock is on the order of 10~5 inches. Due to 

the complex interactions involved, analysis of the pro- 
perties v;ithin the shock are beyond the scope of this course. 

We v;ill analyze the change in fluid properties across the 
shock by investigating the behavior on both sides of a 
shock wave. 

b. Now read Sections ^.1 and H.2 in the text. 

c. We will first consider a "standing (i.e. NON-MOVING) 
normal shock"; that is, a standing shock which is perpendicular 
to flow. With the principles developed in this analysis 

we will then investigate the possibility of the existence 
of an expansion shock. Then moving normal shocks, conver- 
gent-divergent nozzle operation, and oblique shocks will 
be covered. 
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4.2 - Standing Normal Shock 



a. 



Consider a flow situation as shown in Figure *4.1. As 
v;ith previous analyses v;e first establish a control volume. 



FLOW 




Figure *4.1 



then apply the basic equations to the properties across 
the control volume. 



b. The following assumptions will be made : 

1) Since the shock is very, very thin, v;e may consider 
the control volume equally thin and assume that 
there is no surface along the wall to create friction. 

2) Due to the thin control volume, the cross-sectional 
area on both sides of the shock may be considered to 
be essentially equal. 

3) V/e will continue to assume steady, one-dimensional 
flov; with negligible potential changes. 



c. Consider first the continuity equation. Since this is 
steady, one-dimensional flow: 

♦ 

m = constant (*4.2,1) 

or 
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but 



A 



1 



A 



2 



so 



P 






9 



(» 4 . 2 , 2 ) 



d. Now recall the energy equation for steady, one-dimensional 
flow between two points; 



^1-^2 


h, , 

t2 + w 


(*2.1,7) 




®l->2 






We can consider 


this an adiabatic 


process with no 


shaft vjork 


being done, so 


*^1^2 ~ '^s ~ 

®l->2 


0 , or 





h 



t 



1 



9 



(*^. 2 , 3 ) 



e. The x-component of the momentum equation for steady , 
one-dimensional flow (Equation 1.13) becomes; 





(V.dA) 



(*4.2,i|) 



Integrating the right side of equation *4.2,2 yields 




(V.dA) 





We can also express the summation of forces in the x- 
direction as 




PfAi P 2 A 2 



(P^ - P2)A 
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So equation *4.2,4 becomes 






(*4.2,5) 



but m = pAV so show that; 






S. 



?2 + 



P2V2' 



(*4.2,6) 



To this point v;e have been dealing with a general fluid 
By introducing the Perfect Gas Law 



p = 


pRT 


(1.23) 


and recalling the expression for the speed 
a Perfect Gas (Equation *3.3,4): 


of sound in 


V = 


Ma = M yj yg^^RT 


(*4.2,7) 


We can 


now express Equation *4.2,2 in terms 


of Mach number 


P M 
^l"l 




(*4.2,8) 


Vh 


yPh 


Obtain 


this expression yourself for drill. 





Note: .Y is assumed constant 

for this expression. 



Similarly, recalling from Equation (3.6) 



= T[1 + ^ 



(3.6) 



and since 
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from Section 1.1 in the text, the energy equation 
(Equation can be written 






T2CI + ^ 



(*^.2,9) 



g. Similarly, using the definition of Mach Number, the 

Perfect Gas Law and the expression for speed of sound in 
a perfect gas, shov/ that the Momentum Equation (Equation 
(^.2,6) can be written: 



and three governing equations, nos. *A.2,8, *^1.2,9, and 
*^1.2,10. Given any ^ of these we should be able to solve 
for the remaining three variables. 



i. We now find it possible to combine the three basic equations 
so as to derive an expression for M 2 in terms of M^ and 7 . 
Rearrange the continuity and energy'^equations as follows: 



P^[l + yM^^] 



P^Cl + yM^^] 



(*il.2,10) 



h. We have in effect seven variables: 



Y, T^, P^, M 2 and T 2 




(«i|.2,8a) 



1/2 




(^A.2,9a) 



and the momentum equation: 



P 



1 



1 + yM2^ 



(«-4.2,10a) 



P 



2 



1 + yM^^ 



9 ^ 



Now combine these three and shov; that: 



1/2 



~1 + yM2^ 




1 

CM 

CM 

1 — 1 
1 CM 

t — 1 

( 


1 + yM^^ 


, ^2 


1 + 



Solve for and obtain: 



M 



2 



2 




+ 



)M 



2 

Y-1 

2 

1 



1 



(»^.2,lla) 



So for a flow situation such as in Figure *4.1, knowing 
the conditions at (1) (i.e., before the shock) and solving 
for M 2 we can utilize equations *4. 2, 10a; *4.2,7; *4.2,6; 
*3»3,15; *4. 2,9a; and equations 3-6 and 3-7 from the text 
to solve for all other properties. 

I REPEAT ! If the properties are knovm ahead of the ) 
shock (condition 1) then ALL conditions ^ 
after the shock (condition 2) can easily ( 
be found. ( 

For a given fluid (i.e., y is known) since the property 
ratios are functions of y and only , it is possible to 
pre-compute the ratios of the fluid properties, and 
tabularize them. Appendix B in the text is such a table. 
Take to be 2.0 and solve for M 2 , and the property 

ratios, using the relations we have developed. Compare 
your ansvjers with the values in the Normal Shock table. 



Notice that it is possible to use a subsonic value of Mq 
and obtain a MATHEMATICALLY correct answer for M 2 which 
is INVALID. Why? 



HINT: 



Take M^ to be 



577 and solve for 



t2 



tl 



Now what sort of value must As (in this case 
Asi) have to satisfy Equation *3-3515? This 
proves that an "expansion shock" is NOT 
POSSIBLE. 
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k. Examine Appendix B in the text. Note that as M, increases 
M 2 decreases and 



P 



2 



P 



1 



^1 



T 

2 

and Fp— are all > 1.0 
^1 



Therefore a shock is always a compression process, with 
?2 always greater than Notice that Pt2 is always 

less than P^^ (though for values of M]^ near 1.0 only 
slightly) . 

Thus a shock process is IRREVERSIBLE . 



1. Read Section 4.3 in the text and study the examples. 
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^ . 3 - Moving Normal Shocks 



T = 70 



V 



= 400 f ps' ' ^ 




Figure *^1.2 



P = 10 psig 



valve 



open 




a. Now consider a situation as pictured in Figure ^^.2 where 
air is flowing in a pipe at the conditions shown and the 
valve at the end is suddenly closed. This will cause a 
shock wave to propagate back through the duct as pictured 
in *i|.2a. VJe do not know the speed of this shock wave. 

Is this a steady flov7 problem? If you said "Yes", 




and while you're there, re-read your handout entitled 
Comments on Steady Flow. 
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b. We can change this problem into a steady flow problem by 
"hopping aboard the shock wave." 







V, = 400 + V L / 


Q) V„ = V 


1 s ?N 






r ^ 


= 10 psig y 


p = 9 

2 


T^ = 530 °R 1 


T =9 
^2 



Figure ^k.2h 



We now have a steady flow standing normal shock problem. 

The velocity upstream of the shock wave will be 

and the speed doivnstream v/ill be . Why can we 

not derive a so-called closed form solution for M 2 ? 

c. A "trial and error" process is required to solve this 
problem. Let us assume that = I .3 

a^ = ^ (1.^^) (32.2) (53.3) (530) = 1128 ft/sec 

Thus V^' = M^a^ = ( 1 . 3 ) (1128) = li^67 ft/sec 
and Vs = - 400 = IO 67 ft/sec 



^2 

From the shock tables obtain M 2 = .786 and jj— = I .191 

Thus T 2 = 1.191 (530) = 632 °R 

a 2 = 1233 ft/sec 

and V 2 = M 2 U 2 = 970 ft/sec 
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But V 2 must equal Vg , which it does not! Thus, our 
assumed Mach Number v/as in error. A few trials will 
converge on an assumption of Mx = 1.235* Follow these 
calculations through and show that V 2 = Vs = 992 ft /sec. 



What v;ould the pressure be behind the shock wave? How 
do these values compare with those in the actual "unsteady" 
problem posed in Figure *4. 2a? 



d. Read Sections 4.4, 4.5 and 4.6 in the text. 
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^ ^ - Convergent-Divergent Nozzle Operation 



a. We recall that a convergent section acting as a nozzle 
will become "choked" at M = 1. 



If the pressure downstream of such a nozzle is 
decreased the flow will accelerate until M = 1 
is reached at the exit (see Figure *4.3). 





Figure *4.3 



Further decrease in pressure outside the nozzle 
exit will not induce a further decrease in 
pressure Inside the nozzle. The nozzle is said 
to be "choked". 



We also recall (from Figure *3.1) that in order for a 
device to function as a 'nozzle' above M = 1.0 (l.e., 
continue to accelerate the flow, thus converting more 
enthalpy into kinetic energy) the cross-sectional area 
must begin to increase . Thus we can presage the existence 
of a convergent-divergent section (\'7lth M = 1.0 at the 
throat) which functions as a nozzle throughout its entire 
length. This is called a DeLaval nozzle. 
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De Laval Ilozzle 




Entrance Throat Exit 



First 

Critical 



Third 

Critical 



M=1 



Figure *4.^1 



The convergent section of the C-D nozzle offers little 
design difficulty, but the divergent section must be 
carefully designed to reduce losses due to separation, 
wave interaction and other fairly complex phenomena. 



b. Belov/ the diagram of the C-D nozzle in Figure *4.4 is 

depicted a plot of static pressure ratio versus location 
in the nozzle. V/e can envision tv;o possible flow condi- 
tions with M = 1 at the throat. Each of these conditions 
is associated v;ith a pressure value at the exit known as 
a "critical point." The "First Critical Point" represents 
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flov; which reaches M = 1.0 at the throat, but is subsonic 
in both the convergent and the divergent sections (dotted 
line). The "Third Critical Point" represents the design 
point (or pressure ratio) of the C-D nozzle. Both 1st 
and 3rd critical points represent isentropic flow condi- 
tions. Any exit pressure above 1st critical v/ill result 
in subsonic flov; throughout the nozzle, or typical 
"venturi" operation. A point between 1st and 3rd critical 
will Induce non -lsentropic flow through the nozzle. 



c. Since both 1st and 3rd critical represent isentropic flow. 
Appendix A, "Isentropic Flov; Tables", may be used to 
determine the values of Mach Number at first and third 
critical, provided v;e know the area ratio betv;een the exit 
and the throat . 

As an example suppose: 



^ = 3.0 

^throat 



We know that at the throat M = 1.0 so ~ > or: 

^ = 3.0 

Prom Appendix A in the text : 



for 


first 


critical 




P 

Pt 


« .973 


(*4.4,2) 


for 


third 


critical 




P 

P. 


= .0472 


(*4.4,3) 



I Two Notes on this solution 

1) Linear interpolation can be used to obtain 
better accuracy. 

2) The small figure (-1) indicates the pov;er 
of ten with which the tabular value must be 
multiplied . 

It is also interesting to compare these results to that 
of a converging (only) nozzle. In the case of a converging 
nozzle with y = 1.4 we need a pressure ratio of .528. 
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(Pgxit ^Inlet^ produce sonic velocity. In the above 

converging-diverging nozzle we have sonic velocity (in the 
throat) v;ith a pressure ratio of only .97. 



d. As you might expect, somewhere between "FIRST BASE" and 
"THIRD BASE" there must be a "SECOND BASE". Let's 
consider what happens in the region between "first critical" 
and "third critical." 



e. As the receiver pressure is lowered below "first critical," 
a normal shock forms just downstream of the throat (see 
Figure *4.5). The remainder of the "nozzle" is now acting 
as a diffuser since the flow is subsonic and the area is 
increasing. The shock will locate itself in a position 
such that the pressure changes that occur ahead of the 
shock, across the shock, and downstream of the shock v/ill 
produce a pressure that exactly matches the outlet pressure. 
In other words, the outlet pressure determines the location 
and strength of the shock. As the pressure is further 
lowered, the shock continues to move tov;ard the exit. When 
the shock is located at the exit plane this condition is 
referred to as the 2nd critical point. If the receiver 
pressure is between 2nd and 3rd critical, then a compres - 
sion takes place outside the nozzle. This is called 
"overexpansion" (i.e., the flow has expanded too far v/ithin 
the nozzle). If the receiver pressure is below 3rd critical 
then an expansion takes place outside the nozzle. This 
condition is called "underexpansion." We shall investigate 
these conditions later in the course. 




Figure *4.5 
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f. We can now use what vje know about standing normal shocks 
to find the conditions at '2nd critical'. Just ahead 
(e a\vay", as the mathematicians are vjont to say) of the 
shock at 2nd critical, the conditions are identical to 
those at 3rd critical (or design). If you have not done 
so, use the tables to confirm the exit Mach at 3rd critical 
is M = 2.6^ for an Area ratio of 



A 

•^exlt 

^throat 



3.0 




Shock — 

2nd critical 



Figure *^.6 



P . . Po 

exit _ 2 

P P 

^inlet tl 



If we assume stagnation 
conditions at the inlet. 



Is this a reasonable assumption? 
From the Tables : 




tl 

L 



From Isentroplc F'low Tables § M = 2.6^ 



From Normal Shock Tables @ M = 2.64 
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Compute this and verify that: 



p 

exit 

^inlet 



.376 



g. Now suppose vie subject the same De Laval nozzle to a 

pressure ratio of 0.^, whre will the shock form? Is this 
* underexpanded ' or 'overexpanded' Plow? Think about it 
before checking the solution below. 

Solution 

Since the flov; is between 1st and 2nd critical, a normal 
shock has formed Inside the nozzle — hence the terms 
"underexpanded" and "overexpanded" do not apply to the 
flow . 




inlet 



exit 



We have nov; enough information to analyze the flow, 
including losses. 

Known 



exit 

■^throat 


= 3.0 


exit 


- O.i} 


^inlet 
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But : 


^inlet 


inlet 






and 


■^throat 


A*. ^ ; 

inlet < 


[ since 
1 losses 
1 inlet 


we assume no 
between the 
and the throat. 



Now recall Equation ^3 -5, 3- Therefore: 



''iJlet ” *eJlt h " Constant (‘ll.ll.iO 

exit 

Thus ; 



- 




- — 




— 




- - 




— “ 




” 


A .4. 

exit 




P P.. 
exit 




A 

exit 




P - 4 - 

exit 


= 


A . 4- 
exit 




P • 4- 
exit 


■^throat 




P- T ^ 
inlet 




A ^ 

^inlet 




_ ^inlet_ 




A * 
exit 




L exitj 



t f 



knovm known 

I I 

(3.0) (O.iJ) = 1.2 



or: 



~ 




— 


A . 4 . 
exit 




P -4- 

exit 


A * . , 




P 


exit 




^exit 



Prom the equations developed in Unit 3 ’ 



A 



f (M, y) 



(*3.4,2) 



= f(m,y) • (3.7) 

^ t 
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You can obtain an expression for 



A 



_P. 



= f(M,y) 






John does not tabulate this value but a sample table 
including this function was passed out in class. Excerpts 
from a similar table for Y = 1.^ show: 



M 



A -L 

Pt 



.47 


1.205 


.48 


1.178 



Interpolation reveals the Mach No. at the exit to be about 



M 



exit 



. H72 



If we can determine the stagnation pressure ratio 

P, /p v/e will be able to locate the shock, 

^exit' ^inlet 



exit 

^inlet 



P 



'exit 



P 

exit 



exit 



P, 

^Inlet 



1 . 


' 1 


It 


1 


1 .859 / 


• T 



\ 

' Given 

‘From Isentropic tables for M = .472 
Since all losses occur across the shock we can say that 



t . . 
exit 



'inlet 



(«4.4,6) 



across shock 
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From the shock tables find M, = 2.59 and from the 
Isentropic tables find 



^shock 

^throat 



2.869 



We also have a measure of the loss from 



P 



t 



2 





(*3.3,15) 



As we recall that AS = as^ here. 



h. Read Section 5.2 and 5.3 in the text. 
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4.5 - Oblique Shocks 



a. Read Section 6.1 in the text. 

b. Consider a Standing Normal Shock 



Recall: using Shock as 

Reference 

Vin > Sonic 
V^n < Sonic 

Figure ^4.8 






V 



In 



V 



2n 



Now, let us superimpose a velocity to & V2n (say Vt) 

Note that this is the equivalent of running along the 
shock front. Recall that by doing this we have not 
changed the static states of the fluid, but we have 
changed the stagnation conditions. We then have: 



V 



n 





Since ¥ 2 ^ 



V 



In 



and 



Notice that V 2 has turned 
towards the shock front. 



Figure ^4.9 
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We normally picture this as follows: 




Flow 

Deflection 

Angle 



Figure *^.10 



NOTE (1) Flow direction has changed 

(2) > V 2 as with normal shock 

Vj must be supersonic (since is supersonic 

V2 can now be supersonic — if V-j; is large 
enough 

It should be clear that the Normal Shock Relations (and 
Tables) can be used AND ONLY IF , proper care is taken. 



Dividing both sides by a^, and noting that a^ = we 

have : 



V 



In 



a 



In 



— sin 6 
^1 
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Recalling the definition of Mach Number, this becomes 
M^^ = sin 0 

Static pressure ratios, static temperature ratios, etc. — 
can now be taken from the Normal Shock Tables, since these 
were unaltered by the superposition of V, on the original 
shock picture. 



c. Nov;, we know that ^ 1 



thus : sin 0^1 



Let us consider the range of possible 0 for a given M: 
The minimum 0 v;ill occur when sin 0 =1 



0 



min 



M, 



Recall that this is the same expression that was developed 
for the Mach Angle y. Hence, the Mach Angle is the 
minimum possible shock angle. Note that this is a limiting 
condition and really no shock exists for this case since 
then Mij^ = 1.0. The maximum value that 0 can achieve is 
obviously 90°. This is another limiting condition and 
represents our familiar normal shock. 

Notice that as the shock angle 0 decreases from 90° to the 
Mach angle y, M]_j^ decreases from M^ to 1. Since the 
strength of a shock is dependent upon the normal Mach 
number we have the means to produce a shock of any 
strength equal to or less than the normal shock. 



c. Study Example 6.1 on pages 105-106 of the text. 

d. We will nov; try to relate the deflection angle (6) to the 
shock angle (0) starting with the continuity equation. 



^1 In 



= ^2^2n 
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so 



ii 

P2 



V 



2 n 



V 



In 



(*^. 5 , 5 ) 



and from equation (*^1.5,1) and figure 6.5. 

cos 0 = V2 cos (0 - 6) (*4.5,6) 

so 



p -I 

1 cos 0 



cos (0-6 ) 



sin 0 
sln( 0-6 ) 



tan 9 
tan ( 0-6 ) 



(* 4 . 5 , 7 ) 



From our Normal Shock Relations in Unit 4 . 2 , we can show 

that for a Perfect Gas (as a function of M., ): 

In 



^1 (y- 1 ) 2 



(*4.5,8) 



but sin 9 



(*4.5, 8 a) 



Combining equations * 4 . 5 , 7 , * 4 . 5,8 and * 4 . 5 , 8a yields: 



(y+ 1 ) sin^ 6 _ tan 0 

(y- 1 )M^^ sin^O + 2 tan( 0 - 6 ) 



So 0 = f( 6 , M, y) 

or for a given and y, the shock angle (0) is a function 
of the deflection angle (6). But note that this is a 
transcendental function. So for supersonic flow along a 
wall v;e cannot solve for 0 as a function of 6 , M and y 
EXPLICITLY, but we can obtain an explicit solution for 
6 = f(0, M, y). 
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tan 6 



2 cot 0 



1 



(*4.5,10) 



M, 2 . 2^ 

1 sin 0 - 



(y+cos 20) + 2 



Hence we can construct a chart such as in Appendix C to 
relate these three in accordance with Equation *4.5,10. 



e. Now read the remainder of Section 6.2, and carefully 
study the examples. 



f. Read Section 6.3 in the text. 



g. Read Section 6.5. We summarize here some of the 
important concepts : 

1. Flow always turns "toward" an oblique shock front. 

2. For given values of 6 and M, two values of 0 may 
exist . 

a. A large pressure ratio results in a strong 
shock and subsonic M 2 . 

b. A small pressure ratio results in a weak shock 
and supersonic M 2 . 

3 . A maximum value of <S exists for a given Mach Number. 

4. For 6 >6max , a detached shock results (see Figures 
6.8 and 6.9 in text). 
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UNIT 5 - PRANDTL-KEYER FLOW - OBJECTIVES 



The student shall be able to: 

1. Know how Entropy Changes and Pressure Ratios vary v/ith 
deflection angles. 

2. Explain hov; finite turns (with finite pressure ratios) 
can be accomplished isentropically . 

3. Show Prandt 1-Meyer flovj (both expansions and compressions) 
on a T-S diagram. 

Demonstrate the development of the relation betv/een 
MACH number (M) and flov; turning angle (v). 

5. Shov; hov; tables can be developed for Prandt 1-Meyer Flow 
by the introduction of a reference state. 

6. Explain the governing boundary conditions and shov; the 
results when shock v;aves and P-M v;aves: (a) reflect off 
physical boundries, and (b) reflect off "free" 
boundaries . 

7. Draw the wave forms created by flow over rounded and/or 
wedge-shaped v;ings as the angle of attack changes. Be 
able to solve problems of this type. 

8. Describe and sketch what occurs as fluid flovjs past a 
smooth concave corner and a smooth convex corner, 
explaining why an expansion shock CANNOT occur. 

9. Explain v;ith the aid of diagrams the flow conditions at 
the outlet of a supersonic nozzle created by under- 
expansion and overexpansion. 

10. Solve typical problems using PRANDTL-MEYER flow tables. 

11. Demonstrate by a sketch, an understanding of the 
operating conditions of a fixed-geometry diffuser. 



UNIT 5 - PRANDTL-MEYER FLOW 



5.1 - Introduction to Prandtl-Meyer Flow 



a. You saw from the "Special Problem" that for values of 

(in a NORMAL shock process) v;hich are close to 1.0, the"^ 
value of the quantity _ 1 )^ v;hlch we will call 'm' , 

is small. It was shovm that for small values of 'm' , 
the pressure change across the shock (AP) is proportional 
to the first pcvfer of 'm' . We also sav; that the entropy 
change ( As) is proportional to the third power of 'm' . 



b. Now consider equation *4.5,9 from the last Unit, but 
Inverted. 



tan( 9-6 ) 
tan 0 



(Y-l)M^^sin^e + 2 
(y+DM^^sin^ e 



(* 4 . 5 , 9a) 



We can rearrange this to show: 



2 2 
sin e 



(Y+l)tan(e-6) _ ( Y-1) 
2tan 0 2 



(* 4 . 5 , 9b) 



By applying the auspicious trignometric relations and 
re ranging vje can show that this becomes: 



M^^sin^0 - 1 



(y+1) m sin 0 sin 6 
2 cos(0-6 ) 



(* 5 . 1 , 1 ) 



Notice that for small values of 6 where we can approximate 
sin 6 as 6 , and cos(0-(J)) as cos 0: 



M^^sin^0 - 1 



tan 




(* 5 . 1 , 2 ) 



You saw from equation * 4 . 5,8 as applied in the "Special 
Problem" : 




2y 

Y + 1 




1) 



(*5.1,3) 



115 



Note that in equation *5.1,3 represents the normal 
mach number since from unit ^.5(b) vie saw 



sin 6 



Thus Equation (*5.1,3) as applied to oblique shocks is: 



- P. 



(M^^sln^e - 1) 



(*5.1,i0 



Thus we can now see that for small deflection angles (and 
hence weak shocks) combining Equations (*5-l,^) and 
(*5.1,2) yields: 



P 



2 







(*5-1,5) 



c. We can draw the conclusion from this that for very weak 
oblique shocks the pressure change across the shock is 
proportional to the deflection angle. 



summary: 


AP 


m 


where : 


m = (M^^^ - 1) 


But also 


AP 








and 


As 


« m 







Thus 



As « 
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Now suppose we let the deflection angle become infinitely 
small. What happens to As? If as is essentially zero 
for an inf intesim.ally small turn, vie can accomplish a 
finite turn isentroplcally by a series of such infinites- 
imally small turns, and the finite turn will have a 
finite pressure change. 

Furthermore, since this situation has no losses, it 
is reversible . Thus v;e can have compressions or expansions 
depending on the boundary conditions. 
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d. Nov; read Sections 7.1, 7.2 and 7*3 in the text. Note 

carefully the reasoning which leads to the Prandtl-Meyer 
expansion fan pictured in Figure 7. 8... and it is ISENTROPIC. 
In this isentropic (or Prandtl-Meyer) expansion process: 

1) Mach number Increases . 

2) Mach angle Decreases . 
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5.2 - Flow Equations 



a. Read section 7-^ in the text, through page 125. 



b. Notice that the tables are referenced to the boundary 
condition where Mq = 1.0 and v = 0. Also notice that M 
is the Mach number AFTER the turn. 



c. Nov; consider Figure 7.10 in the text. Suppose you were 
given M]_ = 2.13 and asked to find M 2 , given the angle 
V 2 ~vi = ^0°. From the Prandtl-Meyer flow tables in 
Appendix D, you can find the "virtual" angle i , which 
would have produced a Mach number of M^^ = 2.13. If we 
then add this value ( 30 °) to the value of vi-vp, and 
consult the table for v = 70 °, we find M 2 = 5.3^. 



d. Since this is an isentropic process, we know that Ptl=Pt2 
and T^]_=T-t 2 5 so we may utilize the ISENTROPIC tables 
(appendix A) to solve for static pressure and temperature. 



e. Now read the remainder of Section 7.^ in the text and 
follow through the examples carefully. 



f. Read Section 7.5 in the text. Note that for the 

compression process the assumption of Prandtl-Meyer 
flow (hence Isentropic) is only valid in a limited region 
close to the wall. The discussion on concave corners 
(Section 7.3 in the text) shows that for many situations 
involving a smooth concave turn, the flow characteristics 
away from the boundary or wall can be found by considering 
the turn as a sharp concave corner. The transition from 
point to point close to the wall approaches the condition 




Figure ^5*1 
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of a series of Mach waves as shown in Figure 7*5j but as 
we consider points av/ay from the wall, v;e find the 
inf init esmally small compression waves coalesce into a 
finite oblique shock. Figure 7.6 in the text is not an 
accurate d piction of the coalescing process. If the 
area close to the wall is enlarged, as in Figure *5.2. 




A 



We find an envelope of Mach lines forms as the Mach lines 
coalesce forming a finite shock wave at the Inclination 
Angle, (e) corresponding to the overall deflection (6) 
and the Initial Mach number. Since boundary layer and 
other real gas effects at the wall are often neglected 
in basic engineering considerations, v:e often consider 
only the flow av/ay from the wall. Hov/ever if the turn 
is sufficiently smooth and gradual, the Prandtl-Meyer 
flow region may be significant, and the flow within a 
reasonable distance from the boundary may be treated as 
Prandtl-Meyer flow with reasonable accuracy. 
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5.3 - V/ave Reflections 



a. Read Sections 7.6 and 7.7 in the text. 



b. Review the summary in Section 6,5 of the text which puts 
great emphasis on the Importance of boundary conditions. 

Let us now consider the two principle physical conditions 
that govern wave behavior. We shall accomplish this by 
analyzing the operation of a converging-diverging axl - 
symmetric nozzle depicted in Figure *5.3 v^hich is operating 
betv7een its SECOND and THIRD critical points. Obviously 
at a physical wall, the flow must be parallel to the wall. 
V/e may view the central streamline as though it were a 
"viall" , hence (if we treat the flow as one dimensional 
we may consider the flow in region 1 to be parallel to 
the centerline as depicted. For the second condition 
consider the boundary of the "free jet" as depicted by 
the dotted line in region . There must be pressure 
equilibrium along this boundary, so ?2 - Pambient* We 
can now follow from region to region and by matching the 
required pressures or flovj direction angles, v;e may 
determine a great deal more about the flov; conditions. 




Figure *5.3 
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c. If we operate this nozzle with a pressure ratio betv;een 

2 nd and 3 rd critical it is obvious that we need a compres- 
sion process at or beyond the exit in order for the flow 
to end up at the ambient pressure. However, a normal 
shock at the exit v/ill produce too strong a compression. 
What is needed is a shock process that is weaker than a 
normal shock and the oblique shock has been shown to be 
just this. Thus, at the exit v;e observe oblique shock 
(2) at the appropriate angle so that = P 



d. V/e recall that the flov; across an oblique shock is alviays 
deflected toward the shock and thus the flov; in region 
is no longer parallel to the centerline. Wave front @ 
must deflect the flow back to its original axial direcrEion. 
This can easily be accomplished by another oblique shock, 
(an alternate way of viewing this is that the oblique 
shocks from both upper and lower lips of the nozzle "pass 
through each other" v;hen they meet at the centerline. If 
one adopts this philosophy one should realize that the 
waves are altered in the process of traveling through one 
another . ) 



e. Now, since ?2 ~ Pamb> passage of the flow through oblique 

shock @ v;ill make P 3 > Pamb and region can not have 

a free surface in contact v;ith the surrour^ings . Conse- 
quently a wave formation must emanate from the point where 
wave meets the free boundary and the pressure must 

decrease across this wave. An "expansion shock" would 
fulfill our requirements but we know that no such animal 
can exist. We now realize that wave form @ must be a 
Prandt 1-Meyer expansion so that P^^ = 

f. Hov;ever, passage of the flow through the expansion fan @ 

causes it to turn away from the centerline. Thus, as each 
wave of the P-M "expansion fan" meets the centerline a 
wave form must emanate to turn the flow parallel to the 
axis again. If v;ave were a compression in which 

direction would the flow turn? We see that to meet the 
boundary condition of flow direction wave (d) must be 
another P-M expansion. Thus the pressure in region 

is less than ambient. 



g. Can you now reason that to get from to and me^ 

the boundary condition imposed by the free boundar;:^ (E) 
must consist of P-M compression waves. Similarly 
must consist of P-M compression v:aves in order to turn 
the flow' in region (?) to match the direction of the 
"v;all". Now, is Py equal to, greater than, or less ^an 
Parnb? should realize that conditions in region (?) 

are similar to those in region (J) and so the cycle repeats. 



121 




{ 



h. The situation described above (operation between 2nd and 
3rd critical) is referred to as OVEREXPANSION. It is also 
described in section 8.3 in the text. 



i. From this example we may dravj some general conclusions 
about reflections: 

1) Reflections from a Physical (or pseudo-physical 
boundary where the boundary condition concerns 
the flow direction) are the same "family". That 
is, shocks reflect as shocks, compression waves 
reflect as compression v;aves, and expansion waves 
reflect as expansion waves. 

2) Reflections from a "free!" boundary (v;here pressure 
equilization occurs) are of the opposite family, 
i.e., compression vjaves reflect as expansion vjaves, 
and vice versa. 

j. The angle through which the flow turns from to (J) 
is determined bv the nressure change reauired. v.'hereas 



paralleling the "vjall". 

k. Let us now examine an UNDEREXPANDED nozzle. This is 
operation below 3rd critical. 










P' 




Figure *5*^ 



122 




1 

i 



< 



I 

'll 

i 



i 



Note that the flovi leaving this nozzle has a pressure 
greater than ambient and the flov; is parallel to the 
axis. Reflect back to Figure *5.3 shovjlng the overexpanded 
nozzle and you will see that the condition is exactly the 
same as region ( 3 ) . Thus the flov; patterns are the same 
from this point on. This situation is also described in 
Section 8.3 of the text. 



l. We should note that at the jet edge in a real gas a 
SHEAR layer exists as a result of the velocity gradient. 

m. Now read Section 7.8 in the text. 
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5 • ^ - Supersonic Diffuser 



a. We will nov/ examine the design of a fixed geometry 
supersonic diffuser. Prom our studies of supersonic nozzles 
we v/ould assume that a converging-diverging section would do 
the trick - and Indeed It v/111. However, there are some 
practical operating difficulties that must be considered. 

b. Suppose we design the Inlet diffuser for an airplane 
that will fly at Mach 1.86. Prom the Isentroplc tables we 
see that the area ratio corresponding to this Mach No. Is 
1 . 507 . We shall construct the diffuser with an area ratio 
(Inlet to throat) of I. 5 OO. Nov/ let us follow the operation 
of this diffuser as the aircraft takes off and accelerated 
to Its design speed. Note that as the flight speed reaches 
Mq ^ .^3 the diffuser becomes choked. (Check the subsonic 
portion of the Isentroplc tables for the above area ratio.) 
This condition Is shown In Plgure *5. 5a. Nov/ Increase the 
flight speed to Mq = .6. "Spillage" of external diffusion 
occurs as Indicated In Plgure ^5 .5b. As Mq Is Increased to 
1.0 there Is a further decrease In the "capture area" 

(area of the flov/ at the free stream mach number that 
actually enters the diffuser). 

c. As we increase Mq to supersonic speeds a detached 
shock wave forms. See figures *5-5d and *5-5e. Note that 
at the higher flight speeds the shock moves closer to the 
Inlet as less external diffusion is required to reach 

M]_ = .43 at the inlet. Also note that It is necessary to 
fly at slightly greater than Mq = ^4 . 19 in order to have the 
shock attached to the inlet as shown In figure *5-5f 
(Check the shock tables to substantiate this.) If we now 
increase Mq to ^.2 the shock moves very rapidly past the 
throat and forms in the divergent section dov/nstream of 
the throat as shown in figure ^5.5g. This is referred to as 
"swallowing the shock" and the diffuser is said to be 
"started". Under these conditions we no longer have Mach 
1.0 in the throat. (Can you compute the Mach number that 
exists in the throat?) We can now slowly decrease the flight 
speed to the design condition of M = 1.86 and the shock will 
move to a position just downstream of the throat and occur 
at a Mach number of just slightly greater than 1.0. (Thus 
we have a very weak shock and negligible losses.) See 
figure 5.5h. 

d. Two comments can now be made on the above performance: 

1) In order to "start" the diffuser, which was 
designed for 1.86, it is necessary to "over- 
speed" the plane to a Mach number of 4.2. 
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If the plane slows just slightly below 
its design speed — or perhaps minor air 
disturbance might cause it to drop belov'/ 

1.86 — the shock will pop out in front of 
the inlet and the diffuser must be "started" 
all over again. 



M,=.43 

^ 



in (a) thru (f) 

M2= 1.0 



(i) 
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(f ) 




e. It should be obvious that either of these situations 
cannot be tolerated and for this reason one does not see 
fixed geometry converging-diverging diffusers. V/e can 
not summarize the behavior of fixed geometry diffusers 
as is depicted in figure *5-6. 
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5.5 - Airfoils 



a. Read Section 8.5 In the text. 

b. Follov; carefully through the examples In section 8.5* 
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UNIT 6 - FANNO FLOW - OBJECTIVES 



The student shall be able to: 

1. List assumptions and restrictions applicable to FANNO 
FLOW and relate these to real physical situations. 

2. Sketch a Fanno line in the h-s and T-s planes, identifying 
the regions of sonic, supersonic and sub-sonic flow. 

3. Describe the variation of static and stagnation pressure, 
static and stagnation temperature, density and velocity 
as flov; varies along a FANNO line (for both subsonic 

and supersonic flow). 

State what is meant by "choked flow". 

5. Describe the effects of changing duct length in a choked 
Fanno flovj situation in supersonic and subsonic flow. 

6. Define "friction factor", "equivalent (or hydraulic) 
diameter", "absolute and relative roughness", "absolute 
and kinematic viscosity", and "Reynolds number", and 
know how to determine each. 

7. Starting with basic principles, derive expressions for 
property ratios such as T 2 /T]^, P 2 /P 1 J etc. in terms of 
Mach number and specific heats, in a Fanno flow situa- 
tion. (for Honors credit) 

8. Describe (including h-s diagram) hov; the Fanno tables 
are developed by use of a reference location. 

9 . Demonstrate the ability to solve typical Fanno flow 
problems by use of the appropriate tables and/or 
relations . 
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UNIT 6 - FANNO -FLOV/ - STUDY GUIDE 



6.1 - Introduction to Fanno Flov; 



a. We have to this point touched only briefly on the subject 
of friction losses. This unit will analyze the effects of 
flow in a constant area duct with friction. The results of 
this analysis have many direct engineering applications in 
compressible fluid flow, and v;ill pave the v;ay for increased 
understanding of the effects of friction on fluid flov; 
under any circumstances. To simplify the analysis we will 
assume that the flow is adiabatic and no shaft vjork is 
added to, or extracted from the flov.^. This assumption is 
generally quite valid v;hen dealing with reasonably short 
ducts where no special attempt is made to transfer heat or 
v;ork. Additionally, we shall assume steady, one-dimensional 
flow with no appreciable change in potential. This type of 
Flow is known as "Fanno Flovj". 

b. V/e shall first consider the effect of these assumptions 
on the basic equations. Then the general variation of fluid 
properties with Mach number v;lll be considered. Finally, 
working equations for Fanno flov; with a perfect gas will be 
developed, from which we may see hov; values are tabulated 
with respect to a reference condition as has been done with 
shocks and isentroplc flow. Thus v;e vjill simplify the 
business of problem solving. For your ovjn benefit, sumiraarize 
the assumptions that we are making in Fanno Flov;. 
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6.2 - Basic Development of FANNO Flow 



a. Consider first the continuity equation for steady one- 
dimensional flow as developed in Unit 1: 



m 




(* 6 . 2 , 1 ) 



Since the cross-sectional area is constant, this reduces 
to: 



pV 




(* 6 . 2 , 2 ) 



We will assign a new symbol, "G" to the quantity 'pV’ 
which is referred to as the "mass-velocity". So by 
definition : 



pV E G 



(*6.2,3) 



What are the units of "G"? 



b. We will now resurrect the energy equation, which for 
steady flow between any two points in the flow may be 
written in terms of stagnation enthalpy: 



' hj * “1-2 



(*2.1,7) 



but since 



fll-2 ' “ 1-2 ' ° 



then 



’’ti ' 



(*6.2,iO 



We can write this equation in terms of static conditions 
as : 

,2 



V" 



^ 2g, 



gz 






(*2.1,iO 
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Nov 7 assuming the potential remains essentially constant 
and utilizing equation *6.2,3 we can rev;rite equation 
*2.1,4 as: 



h 



+ 



_g1_ 




(*6.2,5) 



Since for a given flow, G and are constant, \ie see that 
for Panno flow h and p are uniquely related. V/e can now 
plot this relationship for various values of G (these are 
called Panno Lines) in the h-v plane (recalling that 
V = 1/p ) : 



S3 > S 2> S I 




Pigure *6.1 



If the fluid is knov/n, one can also plot, as in Pigure 
*6.1, lines of constant entropy (dotted lines) which come 
from thermodynamic considerations. V/e can then replot 
the "Panno lines" in the more familiar h-s plane as shovm 
in Pigure *6.2. 
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Figure *6.2 



c. We may nov; state in general that the combination of the 
continuity equation with the energy equation and the 
thermodynamic equation of state (vjhich established the 
lines of constant entropy) form the locus of points on a 
given Panno Line. 

One factor stands out: 



Since this is ADIABATIC, i.e., dSe = 0, the only 
viay entropy can be generated is by LOSSES, so a 
process can only move toward INCREASING values of S. 
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6.3 - FANNO FLOVJ versus MACH NUMBER 



a. We shall note that due to the irreversible nature of 

FANNO Flow, a FANNO LINE consists of tv;o distinct regions 
separated by a Limiting Point. 




s 



We note that in the upper branch a process (which must 
always result in INCREASING ENTROPY) results in DECREASED 
PRESSURE. Prom Plgui’e *6.1, ve say what as enthalpy 
decreases, so must density decrease, hence velocity must 
INCREASE! 



b. Now consider for yourself what happens on the LOVJER 

branch. Draw an arrow on the Lower Branch of Figure *6.3 
to Indicate which way a process must proceed. As v;e move 
along this branch of the Fanno line: 



Enthalpy 

Density 

Velocity 

Pressure 
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Let us nov; consider the Limiting Point . We are seeking 
to explore a condition where ds = 0. Recalling the energy 
equation (neglecting Potential changes). 



c , 



ht = h H- ^ 


(*2.1,4) 


Differentiate this to obtain: 




0 = dh + 


(*6.3,1) 


®C 


Also recall that 




pv - (j; 


(*6.2,2) 


Differentiating this we obtain: 




pdV + Vdp = 0 


(*6.3, 2a) 


or 




dV = - ^ 


(*6.3,2b) 


p 




Substitute equation *6. 3,2b into Equation *'6.3,1 
for dh: 


and solve 


V*dp 
dh = ^ 


(*6.3,3) 


Sc 




Recall the thermodynamic property relation 




Tds = dh - vdp 


(*6.3,4a) 


Tds = dh - ^ 


(*6.3, 4b) 



13 ^ 



Combining Equations *6.3j3 and *6.3j^b v/e find that 



Tds 



V^dp ^ 

gcP ■ P 



(*6,3,5) 



This expression is valid at ANY point on a Fanno Line. 

We now apply equation *6,3,5 to the limiting point where 
ds = 0, then we find at the limiting point: 

^ = V^dp 

P " g^P 



or 



V 



2 



g 



c dp 



At 

Limit 

Point 



g 



IR 

c 8p 



s=constant 



2 

We recognize this as the expression for a (see section 
2.3 in the text). Hence at the limiting point: 



V = a 



or 



M = 1.0 
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6.k - WORKING RELATIONS for PANNO FLOW In a PERFECT GAS 



a . 



Prom the continuity equation in steady, one-dimensional 
flov; in a constant area duct; 



G = pV = (^ 



(* 6 . 2 , 2 ) 



Prom the Perfect Gas Lav;: 



P 



_P_ 

RT 



and from the definition of Mach number: 



V = Ma 

and for a Perfect Gas 

a = 



Equation (•^6.2,2) can be re-written: 



PM 








b. Now we turn to the Energy Equation. \-Ie showed that 
between tv;o points in a Fanno Line: 



h 







(*6.i^,2) 



Recall that for a perfect gas; 
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So 



T 



t 



1 









We also saw that for a Perfect Gas: 



= T (1 + ^ M^) 



(3.6) in text 



so between two sections 



Tfd + ^ 



T2(1 + ^ ^ 2 ^ 



(* 6 .d 5 a) 



or as a temperature ratio 




1 + 



1 + 



Y-1 

2 

Y-1 

2 




(* 6 d, 5 b) 



c. By expressing the relation in Equation *6d,l betv/een two 
points in the flow we can write 



PlMi 



PpMp 






(*6d,6a) 



This can be restated as a ratio of the static pressures : 




(*6.i|,6b) 
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or, utilizing equation *6.i|,5b we can now obtain an 
expression for a ratio of the static pressures in terms 
of y and Mach Number: 





1 + 



L 1 + 



Y-1 

2 

Y-1 

2 




1/2 



(*6.1|,7) 



Have you seen some of these equations before? For vjhat 
kind of flow? What assumptions have you fed into these 
equations so far? 



d. Now recall the property relation 



Tds 



dh 



P 



What restrictions, if any, are attached to this expression? 
If we now recall from Section 1.11 in the text that: 



dh = c dT (*6.i|,8) 

and 

= R [ -\ ] (*6.i|,9) 

P Y-1 



we can show that 



Tds 




] RdT 



P 



(*6.^,10a) 



and we can rearrange this to show 



ds _ r _x_ 1 ^ ^ 

“R* *- ^-1 T " P 



(*6.^,10b) 
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Nov; integrate Equation *6. 4, 10b between two points in the 
flow to shov; that 



S 




1 




(^6. 3, 10c) 



Utilizing the expressions we have just developed for the 
static pressure ratio and the static temperature ratio 
we can nov; express the entropy change as a function of y 
and Mach Number : 





1 

C\J 

H- 

1 


- In 


1 + ^ 


H 

t — 1 

1 

1 


1 + 


1 + 




2 2 




2 2 



1/2 



(*6. 4,11a) 



You can nov; combine the terms of this expression to show 
that 



S^-S 

R 



1 




1 + 



1 + 



xzi 

2 

Y-1 

2 




Y 

Y-1 



1 

2 



(*6. 4,11b) 



or 



S2-S 



1 




1 + 



1 + 



Y-1 

2 

Y-1 

2 




Y+1 

217^ 



(6.4,11c) 



Refer back to the section where v;e developed the stagnation 
pressure ratio as: 



P 



t 



2 



P 



t 



1 




(*3.3,15) 



You should satisfy yourself that this relation applies to 
Panno Flovi;. 
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If we now convert Equation to; 





1 + 



1 + 



y-1 

2 

Y-I 

2 




Y + 1 

2 ( y-1) 



(*6. ^ ,12) 



Now take the anti-log of both sides to obtain 



Y + 1 

(*6.il,13) 



VJe now have the means to obtain all the properties at a 
downstream point @ if v;e know all the properties at 
some upstream point 1 and the Mach Number at point @ 
But we need to be able ^ predict the Mach Number at 
from the physical set-up of the flow (i.e., Duct lengi^, 
etc . ) . 



R ^2 



M, 



1 + ^ 

1 + - M ^ 

11- 2 “i 



e. In Unit 1, we developed an expression for the Momentum 

Equation as it applies to steady, one-dimensional flow in 
any fluid 



p 



+ 



fdxV^ 



dz 

Sc 




(*1.6,l6) 



If we now apply the perfect gas law and the definition 
of Mach Number to this equation, we can show (neglecting 
Potential change); 



^ (RT) + 



fdx M^ygpRT 

D 2g„ 

e c 



MVygpRT' dV 
Sc 



1^40 



0 



,liJ) 



1 

i 



1 









This can be rewritten as: 



dP fdx M^y 

P 2 

e 



+ ^ dM^ 




2 ^ 
T 



0 



(*6.4,15) 



These two terms result from 
eliminating dV in the last term 
of Equation *6.4,14. Can you 
show this? 



f. VJe must first evaluate dP/P as a function of Mach Number. 
Prom continuity: 




(*6.4,6) 



Therefore : 

In P + In M - ^ In T = In (|; 



(*6. 4,16a) 



Differentiate to shov/: 



dP + ^ _ 1 ^ 

P M 2 T 



( *6 . 4 ,l6b ) 



Now we must obtain an expression for dT/T as a function 
of Mach number. From the energy equation we obtained: 



T^ = T(i + ^ M^) = cj: 

or 

In T + In (1 + M^) = In^ 



(*6. 4,17a) 
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Again differentiate to obtain 



T 



+ 



d(l + 
1 + 




0 



(*6.ij,17b) 



Now combine equation *6.^,17b with equation *6.4,l6b 
and substitute into the momentum equation in the form 
of equation ^'6. ^,15 to show: 



fdx 



D 

e 



( Hi) 

1 + 



2dM 

2 

ym 





(*6.4,18) 



We can now integrate equation *6.4,18 from point (l) to 
point in the flow (see Figure *6.4) to obtain: 



f (X2“X^) 

dI 



Y+1 

2y 



In 



1 + ^ 
1 + 



2 2 
,M2 j 



+ 1 ^^9 

^ in -i 
M ‘ 



T 







Figure *6. 4a 



Figure *6. 4b 
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Knowing the duct length (x 2 - xi), flow rate, fluid and 
duct material (hence the friction factor), configuration 
(to obtain De ) and all conditions at one point, we may 
now find the conditions at the other point. 



g. The solution of problems (as in previous analyses) is 

simplified by the introduction of a reference condition. 
The obvious choice is the limiting condition where 
M = 1.0. If vie let point "(T) be the condition (i.e., 

M = 1.0; reached by Fanno fTow ) then: 



T^ = T» 



Point ( 2 ) may be any other point in the flow. Equation 
*6. 4,5b may now be reduced to: 



Y+1 



T 



2 



(*6.4,20) 






1 + ^ 



Using a similar approach, obtain the expressions for; 



P 



9 



(*6.4,21) 



pTT 



and 




9 



(*6.4,22) 
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You should also be able to corroborate that 



p 



Y1 

V 



1 

M 



1/2 

2(1 + ^ 

Y + 1 



(*6.i|,23) 



Now check your expressions for P/P* and P^/P^* by 
solving for any arbitrary value of Mach (subsonic or 
supersonic) and using a value of y = 1.^. Then compare 
your solutions against the tabularized values in Appendix 
E in the text. 



( 




i 



6.5 - Applications of Fanno Flov/ 



a. Applying the reference condition to equation 
(l.e., = 1, = M) yields: 



f(x-x*) 

D_ 




~1 + ^ 6 ^ 

X+I 

2 



^ In m2 




(•6.H, 



Since X* Is always greater than x (see Figure ^6.^b), 

It makes sense to change all signs In Equation ^6.4,24 
to obtain a positive expression which we can simplify to: 



f (x*-x) 
D„ 




r 4^ 1 


. 1 


r 1 


2 


+ — 


- 1 


1 + ^6^ 


Y 





(* 6 . 4 , 25 ) 



The quantity (x*-x) now represents the maximum duct length 
which may be added to the given duct and still maintain 
the same flow. This expression 



f L 



max 



D, 



Is tabulated In Appendix E In the text. 



b. In order to solve problems we must be able to calculate 
friction factor (f) and equivalent (or "hydraulic") 
diameter (De). From Unit 1 you may recall the definition 
of these two quantities. 



f 




D 



e 



P 

circular duct, Is simply the diameter (D). 



For a 



i 

I 






c. The friction factor may seem a bit Intractible since the 
shear force (t^.) Is difficult to measure. However, you 
may recall from your previous courses in Fluid Mechanics 
that friction factor is a function of Reynolds number (Np) 
and "relative roughness" which v;e can define as c/Dg v;here 
"c" is an average value for the absolute roughness of a 
given material. Table ^6.1 gives some typical values of 
"e" (which has the dimension of "feet") for various 
materials. Thus the relative roughness (e/Dg) is a 
dimensionless parameter. 



TABLE *6.1 

TYPICAL ROUGHNESS VALUES 
Material 

glass, brass, copper, lead 

steel, wrought iron 

asphalted cast iron 

galvanized iron 

cast iron 

wood staves 

concrete 

riveted steel 



"e" (in ' feet ' ) 
smooth < .0001 
.00015 
.0004 
.0005 
.00085 
.002 
.01 
.03 



d. Similarly, Reynolds number 'is a dimensionless parameter 
which relates size and viscosity. Reynolds number is 
defined as; 



Nr 



If we 
u 




where y is "absolute viscosity" 
define a "kinematic viscosity" 

P 
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Then 



Nr 



VD 

V 



Table 3 on page 38? in the text gives typical values for 
for various gases. With this information, and 
knowing the conditions at a given point in the flow, we 
can calculate the density and velocity (remember G = pV 
= (^ ) hence Reynolds number. Then using experimentally 
derived curves such as depicted in Figure *6.5, the 
friction factor may be determined. 




e. Now read sections 9.1 and 9*2 in the text. Study care- 
fully example 9.1. 



Now consider the following example problem: 

Given a circular duct with these specifications 

100 ' 



Diameter - 6 inches 
Length - 100 feet (Ax) 
Material - Galvanized Iron 



(D 



h 



max 



max 



lil7 



Figure *6.6 



It is desired to have this duct deliver 800 cubic 
feet/minute of air at 70®P and l8 psia. Find the inlet 
conditions required at point . 

First, let us suimnarize a step-by-step problem solution 
technique : 

1) Sketch the physical problem. 

2) Establish and designate locations where conditions 
are known and/or desired. 

3) Compute the Equivalent Diameter. 

4) Find friction factor. 

5) From equation ■^'6.4,19 (or the Tables using FL /D, 
if known) find Mach number(s). 

6) Use relations or Tables to find additional properties 
desired . 



Now let us apply this approach to our example problem. 
Steps 1) and 2) are done in Figure *6.6. 

Recall for a circular duct that the Equivalent Diameter 
is the pipe diameter. 



Dg = .5 feet 

From the Table *6.1 for Galvanized Iron the value for 
e/D is 



c_ 

D 



9 



To find Reynolds number v;e must know p and V. Since we 
know the conditions at , we can com*pute 




^2 

P2 



A 



2 



Tj = Il9^f = 1128 ft/sec 

P 

^ = .0918 Ibm/ft^ 

A ip 

ttD^ tt « , 2 

T“ - IF 



li|8 



^ 2^2 ~ Volumetric flow rate 

So 

Vo = ft/sec 

^ ? 

^2 

Mo = -^ = .602 ^ .60 

^ 3. 2 



e 



Now you should be able to determine from the attached 
"Friction Factor for Pipes" chart; 



f = .0197 ~ .02 



So 



fAx _ (.02)(100) 

D 



-From the 


Appendix 


E, for Mo = . 6 , max 

D 


= .^49081 


thus 


f L, 


f Lp 

= + 


fAx 




D 


D 


^ ? 





Now consult Appendix E to find M, = 

? 







i 



Now knowing M]_ and M 2 you can use the property ratios 
to obtain the properties at point (l) • 



T 



1 








556 °R 



9 



Note that vje do not have to bother v^lth subscripts on 
the * quantities (such as and P*) since for any given 
Panno Line there is only one reference location. 



g. Now read section 9-3 in the text. Pay close attention 
to the flovj situations as described in Figure 9 .16. 

We will not delve into the subject of Iso-Thermal Flow, 
but if time permits, you are encouraged to finish reading 
Chapter 9 in the text. 
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Unit 7 . 



RAYLEIGH FLOW - OBJECTIVES 



The student shall be able to: 

1. State the assumptions and restrictions utilized in 
the analysis of constant area flow vjlth simple heat 
transfer (Rayleigh Flow). 

2. Given the general equations of energy, continuity, 
momentum and state, demonstrate hov-; the assumptions in 
objective (1) modify or eliminate each term. 

3. Sketch a Rayleigh line in the P-v plane together v;ith 
lines of constant entropy and lines of constant temper- 
ature. Indicate the points of maximum temperature 

and maximum entropy, and the directions of increasing 
entropy and temperature. 

4. Sketch a Rayleigh line on a T-s diagram, indicating 
the regions of Sonic, Subsonic and Supersonic flovj. 

5. Describe the changes in fluid properties which occur 

as one moves along a Rayleigh line for the case of heating 
and the case of cooling in both subsonic and super- 
sonic flow. 

6. Correlate the reference points of maximum temperature 
and maximum entropy on a T-s diagram of a Rayleigh 
line with those on a P-v diagram. 

7. Explain what is meant by "thermal choking". 

8. Explain by a T-s diagram ho\>j the flow adjusts to the 
addition of heat in a constant area duct in a "thermally 
choked" flow situation in both subsonic and supersonic 
flow. 

9. Solve typical Rayleigh flow problems by the use of 
tables and equations. 
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Unit 7 . 



RAYLEIGH FLOW - STUDY GUIDE 



7 . 1 Introduction to Rayleigh Flov'j 

a) Constant area flov; with significant heat transfer, 
addition or removal, is called "Rayleigh flow". We 
will approach the analysis of this flov; condition 
in a manner similar to that used in previous flows 
such as normal shocks, Fanno floxv, etc. We vjill first 
develop the basic equations as applied to the Rayleigh 
flow situation. Then perfect gas relations will be 
introduced. Finally, we will introduce a reference 
state which permits tables to be developed, which will in 
turn simplify problem solution techniques. 

b) Read Section 10.1 in the text. 
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7.2 - Basic Relations 



a. We can suminarize the assumptions made in Rayleigh flow 
as follovjs : 

1) Steady flov; 

2) One-dimensional flovj 

3) Negligible change in potential energy 

4 ) No work (shaft) 

5) Constant area 

6) Negligible friction 

This last assumption v;arrants some comment since internal 
irreversible effects (dSi) are obviously present. What 
is really being said is that since we are dealing with 
significant heat transfer the change in entropy from this 
cause (dSe ) will be relatively large. Thus the assumption 
that is actually being made is that: 



dSg >> dS^ 

and we can thus neglect dS^ and say that 



dSg dS (*7.-2,l) 

b. From the continuity equation, we v;ill recall 
pAV = m = (f 
but A = (^ hence we can say 
pV = (j: 



Do you recall this expression? Do we have a name for 
this constant? 



c. The energy equation for steady, one-dimensional flow can 
be v;ritten as: 



+ q 



l->2 




+ w 



s 



l->2 



(* 7 . 2 , 2 ) 
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For our flow 



0 , and 



s 



l->2 



+ q 



l->2 




(*7.2,3) 



A Note of Caution: 

This is the first major flow category for 
which 'h-t' has not been constant. Consider 
a perfect gas. Will ’T^’ be constant? 



d. Recall our basic momentum equation in differential form: 



W + f ^ ^ + VdV ^ 

^ ^e ^®c ®c ®c 



= 0 



(*7.2,4) 



Under our assumptions this becomes: 



dP + W = 0 

P ge 



(*7.2,5a) 



or 



dP + 




(*7.2,5b) 



From continuity we know that pV = G = constant. Thus 
the momentum equation can be written as 



dP + 



GdV 



= 0 



(* 7 . 2 , 6 ) 



This can easily be integrated to give: 



P + 



GV 

Sc 



0 



(*7.2,7) 
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Or, by substituting for the velocity, this can be 



written as 






P + 


0= 

ScP 


Constant 


(*7.2,8) 


or 

P + 


g2 ^ 
Sc "" 


Constant 


(^7.2,9) 



For a given flov; rate this indicates a linear relationship 
between P and v. Try sketching equation (*7.2,9) in the 
P-v plane. This line is called a "Rayleigh line." 



P 



Hint: 

y = mx + b 



Figure 



*7.1 



V 



e . 



As an aside, rewrite equation (*7.2,7) as: 



PA + OAVLV 



Constant 



or 



PA 



+ 



m V 
Sc 



Constant 



Compare this expression with equation (4.2,5) in the unit 
on shocks. This constant is sometimes called the "Impulse 
function" and is given the symbol "F" or "I" by various 
authors . 

f. If we add lines of constant temperature to our plot 
(this is easy if vie are dealing with a perfect gas) 
we will have a diagram as shown in Figure *7-2. 
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A 







p 




Figure *7-2 



V 



Consider a heating process which moves from A to B along 
the Rayleigh line. If we add more heat, v;e move farther 
along and the temperature increases more. When point C 
is reached, the fluid has reached a maximum temperature. 
Can we move farther along this Rayleigh line? \7hat is 
the limiting factor? 

Recall that as we add heat we are Increasing the entropy 
of the fluid (dSe = dq/T). Thus we had better investigate 
lines of constant entropy. When these lines are plotted 
on Figure *7-2 they are very similar to the lines of 
constant temperature. To get an accurate picture of v^hat 
this plot v/111 look like we must investigate the relative 
slopes of 'T = Constant' and 'S = Constant' lines in the 
P-v diagram. It will be much easier to carry out this 
investigation if v;e consider perfect gases. 

For a 'T = Constant' line 

Pv = RT = constant 

Pdv +vdP = 0 



and 



dv 



P 



(* 7 . 2 , 10 ) 



V 
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For an ' S = Constant' line 
Pv"^ = constant 

v^dP + Pyv^~^dv = 0 



and 



dv 



Py 



,Y-1 



P 



^ V 



(* 7 . 2 , 11 ) 



Comparing equations (* 7 . 2 , 10 ) and (* 7 . 2 , 11 ) and noting 
that Y is always greater than 1.0, v^e see that the 
'S = constant' line has a greater negative slope and 
thus these lines xvill plot as shown in Figure *7.3. 




Figure *7.3 



g. V.'e now see that we can add more heat and go past the 
point of maximum temperature until vie reach the point 
of maximum entropy. Let us pause for a moment and see 
if this appears reasonable. 

Normally vie think of the addition of heat as causing the 
fluid density to decrease. This requires the velocity 
to increase since pV = constant. The increase in velocity 
requires a certain increase in kinetic energy. Thus, 
some of the heat v;hich is added goes into increasing the 
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kinetic energy of the fluid and some increases the 
enthalpy (temperature) of the fluid. As v;e add more heat 
the required velocity increase (and corresponding kinetic 
energy increase) becomes larger and larger. Eventually 
we reach a point where all of the heat energy added is 
needed for the kinetic energy increase. At this point we 
reach a maximum enthalpy (temperature). Further addition 
of heat causes the kinetic energy to increase by an amount 
greater than the heat energy being added. Thus from this 
point on the enthalpy (temperature) decreases to provide 
the proper energy balance. 



h. Can we go past the point of maximum entropy? Not by heat 
addition. An alternate v/ay of stating this is that no 
more heat can be added to the system ^ ^ wish to 
maintain the same flow rate (l.e., stay on the same 
Rayleigh lineTT Thus, we see that a limit point has been 
reached and we say that the flow is "choked." \le next 
proceed to investigate this limit point. 

At the point of maximum entropy the slope of the ’s = 
constant' line is the same as that of the Rayleigh line. 
Recall the equation of the Rayleigh line: 




Constant 



Differentiate to obtain 



dv 



-G 



2 




(* 7 . 2 , 12 ) 



Recall the slope of the 's = constant' line: 



dv 





(* 7 . 2 , 11 ) 



We now equate these: 





(*7.2,13a) 



and solve for 



= g^yRT 



(*7.2,13b) 
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You immediately recognize this as sonic velocity and 
thus the 'limiting point' divides the Rayleigh line into 
a subsonic branch and a supersonic branch (similar to 
the situation in Fanno flow). Subsonic flovj is on the 
left and supersonic flow is on the right of the limiting 
point in the P-v diagram. 



Another interesting fact can be noted to exist at the 
limiting point. Recall the entropy relation from 
Section 1.10: 



TdS = dh - — (1.22) 

P 

The differential form of equation (10.2) may be written 

dP = _ (*7.2,li<) 

®c 



The property relation in equation (1.22) can thus be 
expressed as : 



TdS 



dh + 



VdV 

Sc 



(•n.2,15) 



At the 'thermal choke point' where M = 1, dS = 0, so: 



0 



dh + 



VdV 

Sc 



(* 7 . 2 , 16 ) 



Integrate this 

.2 



h + 



V‘ 

2g. 



to obtain: 
Constant 



(* 7 . 2 , 17 ) 



We have already shown that this is the expression for 
'stagnation enthalpy' (neglecting potential). 
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Since at this 'limiting point', = constant, dh^ must 

equal zero . Thus is a maximum at this point. 

k. The above situations become more clear if the Rayleigh 

line is plotted in the h-s plane along with its associated 
stagnation curves. This has been done in Figure 




Figure *7*^ 



Rayleigh Line 



Now correlate this plot with the previous information 
and satisfy yourself that the static and stagnation lines 
are properly drawn. Compare this figure to Figure *7.3- 
Can you add lines of constant pressure to Figure *7-^? 
Note that for a perfect gas this is equivalent to a T-s 
diagram. In v;hich direction does a cooling process move 
along the subsonic branch? How do heating and cooling 
processes move on the supersonic branch? 
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7.3 - V/orklng Relations for Raylegh Flovj - Perfect Gas 



a. We no\-j proceed to develop property relations of a function 
of Mach number and specific heat ratio. These vjill be the 
working equations used for problem solving. We introduce 
the usual perfect gas assumptions as we have done many 
times previously. 



b. Prom the momentum equation (*7.2,7) we can x^rite : 

r2 



P + 



£T 

g. 



= <1 



(*7.3,1) 



but 



= M^a^ = M^Yg^RT 



(*7.3,2) 



and 



^ RT 
Show that 



(*7.3,3) 



P(1 + yM^) = (j) 



(*7.3,^0 



Thus the pressure ratio between the two points is 
P. 1 + yM.^ 

2 = 1 (* 7 . 3 , 5 ) 

^1 1 + y ^2 



How would you obtain a ratio for 




Is this related to the entropy change as we have 
previously always been noting? 
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c. V/e saw that the continuity equation could be expressed as: 



pV = Constant 

Show that this can be written as 



PM 



= (j: 



Nfr 

thus between two points 



(*7.3,6) 



T. 



^ j 



(*7.3,7) 



Combining this with equation (*7.3,5) show that: 

2 



T- 



1 + 
1 + 



2 M ^ 
2 

2 

M, 



(*7.3,8) 



d. In Unit 3 we saw that 



= T [1 + ^6^ M^] 



(3.6) in 

text 



Now you can shoxv that 



T+- 

^2 


1 + 


2 2 
M2 


1 + Y“1 M ^ 

J. T 2 ^'2 


^t 

^1 


1 + y^2^ 




1 + ^ M^^ 



(*7.3,9) 



e. In summary J considering any tvjo points in the flow system, 
if we knov7 the properties at point and one property 

at point (2) we can compute all other conditions at 
point . 
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f. Can we predict what will occur downstream? 

Yes, If we have Information concerning the heat transfer. 

Prom the energy equation we know that 




Thus In 
and for 

^ 1-^2 



^1-2 = ^t. 



general 

perfect gases: 



= c (T - T ) 

p ^2 




(*7.2,3) 



(*7.3,10) 



If we know the rate of heat transfer we can easily solve 
for T |-2 and then use equation (’*7.3,9) to obtain M 2 ; and 
all the other properties are easily computed. 




g. The above developments assume Cp = constant. 

In some cases where heat transfer rates are extremely 
high with large temperature changes resulting, Cp may 
vary enough to warrant using an average value for Cp. 
If, In addition, significant variations In y occur. It 
will be necessary to return to the basic equations and 
derive new working relations, treating y as a variable. 
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7 . ^ - Reference State and Tables 



a. We again Introduce a reference where the Mach number Is 
unity — as reached by the process under consideration — 
i.e., by Rayleigh flow. We designate this the 
condition. Note, in Figure *7-5, hov; this point is 
introduced in both the T-s diagram and the physical 
diagram, for a subsonic, heating example. 




Figure *7-5 



Complete the follov.'lng diagrams (Figure *7-6) for a 
supersonic, cooling example. 




Figure *7*6 



b. 



In the 


previous equations. 


let 


point 


be the reference 


(or ' * ' 


) condition and let 


(D 


be any 


other point in the 


flov/. 


Then . . . 






Ml = 


1.0 




M 2 = M 




It 


p» ^ Ti = T" 




P 2 = P 


, T2 = T 


Pti = 


Pt‘> Ttl= Tt* 




Ptt= Pt 


> 'Pt2" "^t 



I6^i 



So 



P 



2 




P 



(* 7 .^, 1 ) 



1 



Derive the expression for T/T* : 




9 



(* 7 .^, 2 ) 



In like manner, expressions for the other ratios can be 
obtained and all of these are found tabulated in Appendix 
P in the text. 



c. Pick a value for 'M' and let y be 1.^. Compute the 

corresponding value of T/T* and check your value in the 
tables , 



d. Now read Section 10.2 and study the examples carefully. 
We V7ill cover the remainder of Section 10 in the next 
unit which deals v;ith the correlation of Rayleigh and 
Panno flov;. 



e. Note that in Pigure 10.5 in the text, the mass flow rate 
adjusts to a ne\i Rayleigh line if heat is added to a duct 
in V7hich the flow is "choked." 



f. It will prove profitable to prepare a summary of the 
property changes that are taking place in the subsonic 
and supersonic regimes for the heating and cooling 
processes. Try to fill in the spaces in Table *7.1 
without reference to the tables in Appendix P. Use basic 
equations, special equations and/or diagrams that have 
been developed for Rayleigh flovi. Por each property, 
indicate whether that property is increasing, decreasing 
or remaining constant for each regime. 
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TABLE *7-1 



Property 


Heating 


Cooling 


M < 1 


M > 1 


M < 1 


M > 1 


Stagnation 
Enthalpy 
(and temperature 
if a gas) 










Entropy 










Static 

Pressure 










Velocity 










Mach Number 










Static 
Enthalpy 
(and temperature 
if a gas) 










Stagnation 

Pressure 


1 
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Unit 8 - CORRELATION OP FLOWS - OBJECTIVES 



The student shall be able to; 

1. Compare similarities and differences among Fanno 
Flow, Rayleigh Flov; and normal shocks. 

2. Sketch on the same h-s diagram a Rayleigh Line, a 
Fanno Line and a normal shock process (all for the 
same flow per unit area). 

3. Solve typical problems involving various sequential 
combinations of Fanno Flow, Rayleigh Flow and/or 
normal shocks. 
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Unit 8 - CORRELATION OF FLOWS - STUDY GUIDE 



8 . 1 Introduction . 

a) In the previous developments, we have utilized the 
basic equations relating to fluid flow and have de- 
veloped relations which have considered the effects 
of only one of these 'factors' at a time. We can 
gain a great deal more flexibility in problem solving 
by considering various sequential combinations of 
Fanno Flovr, Rayleigh Flow and Normal Shocks. 

b) Another type of problem involves simultaneous addition 
(or removal) of heat vjlth friction effects of the 
same order of magnitude, and may even also include 
changes in area. Exact treatments for problems of this 
type usually involve tedious numerical Integration 

as explicit solutions can not be obtained. Your text 
touches on some of these problems in sections 
9.5 and 10.5- You are encouraged to read these sections 
after you have completed this unit. 

c) In this unit we v^lsh to bring out some of the similar- 
ities and differences among Fanno Flov;, Rayleigh Flow 
and Normal Shocks. Thus we shall assume steady, one- 
dimensional flow in a constant area duct with negligi- 
ble changes in potential and no shaft vjork. 
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8.2 - Panno Flow v/ith a Normal Shock 



a) Recall from Unit 7 the definition of the "Impulse" 
or "Thrust" Function: 




F 



(*7.2,7) 



We saw that in a shock the two states (before and 
after the shock) displayed the same flov; per unit area, 
the same stagnation enthalpy and the same impulse 
function. To corroborate this, re-examine equations 

and *i^.2,6. 



b) A Fanno Line represents states with the same flov; per 

unit area (m/A), and the same 'stagnation enthalpy' (h-^). 
As the flow passes along a Fanno Line, it is experiencing 

the effects of ? Because of this the 'Impulse 

Function' does not remain constant. Suppose we nov/ plot 
the 'Impulse Function' versus Mach number for a Fanno 
Line. Such a plot is pictured in Figure *8.1. 



c) Note that for every point on the supersonic branch, 
there is a corresponding point on the subsonic branch 
where the Impulse Function has the same value. It should 
be obvious that points A and B satisfy all criteria 
for not only the Fanno Line but also for the Normal Shock. 



d) Thus we can imagine a shock between points A and B 
in Figure *8.1. V/ould this shock process occur from 
A to B or from B to A ? 



F 





Figure *8.1 



M = 1.0 



M 
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e) Consider a constant area duct such as in Figure ^8. 2a 
which Is experiencing Fanno Flovj. V/e can nov; visualize 
the possibility of a shock In the duct as pictured 
betv;een points and . This situation meets the 

conditions of flow along the supersonic branch of a 
given Fanno Line from point (l) to point and flow 

along the subsonic branch of the same Fanno Line from 
to (5) . We can also plot this situation In the 
h-s plane as shovm In Figure *8. 2b. 



Flow 











M < 1 






^ 1 


(; 


1) 0 


D 


1 ^ 


1) 



Figure *8. 2a 




Figure *8. 2b 
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8.3 - Rayleigh Flow with a Normal Shock 



a) A Rayleigh Line represents states vjhich exhibit the 
same flow per unit area and the same Impulse function 
but do not have a constant stagnation enthalpy, due to 
heat transfer . We can picture. a Rayleigh Line and its 
corresponding Stagnation Enthalpy line in the h-s plane 
as shown in Figure *8.3* 




s 



Notice that for every point on the supersonic branch 
of the Rayleigh Line there is a corresponding point on 
the subsonic branch which has the same stagnation 
enthalpy. (Also note that the converse is not alv/ays 
true). Would greater than, less than or equal 

to not)? 

b) As with the Fanno Flov/ and a Normal Shock, we can now 
surmise the existence of a duct experiencing Rayleigh 
Flow and a Normal Shock forming as pictured in Figure 
*8.ila. The corresponding process is shov.'n in the h-s 
plane in Figure *8.^b. Note that the flov.' from point 
to point (j) and the flov; from points (J) to lie on 

the same Rayleigh Line. 
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Rayleigh 

Flov/ 



Rayleigh 

Plow 



Flow 



© 



M> 1 






L 



© 



(i) 



M< 1 



7 ^ 






© 



Figure ^8. 4a 
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Rayleigh Flow, Fanno Flovj and Normal Shocks 



8A 



a) 



Consider nov; the combination of Rayleigh Flov/, Fanno 
Flow and a shock process all for the same mass flow 
per unit area, as pictured in the h-s plane in figure 

*8.5. 




b) We have seen that points and have the same 

mass flow rate, the same Imipulse function, and the same 
stagnation enthalpy; and thus can lie on a Fanno Line, 
a Rayleigh Line, and at opposite sides of a Normal Shock. 



c) From the units on Rayleigh and Fanno flows we saw 

that when the "choked flov;" condition was reached along 
the subsonic flow branch, the mass flow rate re- 
adjusted to provide for flov; along a new Rayleigh 
Line (or Fanno Line) if more heat (or duct length, 
hence friction) was added to the flow. In the super- 
sonic regim.e the read j ustm.ent can (and usually does) 
take the form of a normal shock. In this case we do 
not necessarily have to move to a new Fanno or Rayleigh 
Line since once subsonic flow has been obtained we can 
add much more heat, or introduce much more friction, 
as the case may be. Many problems can be considered 
by matching these three types of flow processes 
sequentially, knowing the flow conditions. Let us nov; 
consider qualitatively two such possible combinations. 

d) Consider the com.bination in figure *8. 6a. Assume that 
only v;hen heat is being transferred may friction 
effects be considered negligible. 
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Panno 

Plov; 



Rayleigh Panno 

Plow Plow 



Flow 



=C> 



CD 



(b 



/ 



(D 



© 



Pigure *8. 6a 



If this entire flow system is in the subsonic regime 
the h-s diagram would be that shovm in Pigure *8. 6b. 




Pigure *8. 6b 



We may describe the process as follov/s: 

1) Prom (l) to @ the flow is along a Panno Line 
hence 'h-^^' and ’m/A' are constant. 

2) Prom ( 2 ) to 'm/A' remains unchanged as does 

'P' (i.e., ' PA + mV/gp'), but due to the heat 
transfer, 'h^' is increasing. 

3) Prom ( 3 ) to (J) the flow is along a different 
Panno Line even though it has the same 'm/A'. 

The reason is that the stagnation enthalpies for 
the two Panno Lines are different. 



e. Now suppose v;e consider the situation in Pigure *8. 6a 
but with supersonic flov/ coming in at point . Plot 
this situation in the h-s plane in Pigure *8. 6c assuming 
that the flo;v is not choked. 
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Figure 
*8. 6c 
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f) V/e can easily Imagine the amount of heat added from 

to (D exceeding that required to "choke" the flov/, 
in which case a shock could possibly form. Many 
Interesting possibilities could result depending upon 
the boundary conditions imposed on the system. Are the 
original stagnation conditions maintained? If the flow 
is subsonic at the exit, what pressure must exist? 



g) It might prove useful to summarize the effect on fluid 
properties of flow in ducts in the supersonic and 
subsonic regimes for the various types of flow. 



Table *8.1 





Rayleigh Flov/ 


Fanno 


Flow 


Normal 


Property 


Heating 


Coo] 


Ling 






Shock 




M<1 


M>1 


M<1 


M>1 


M<1 


M>1 


M>1 


Pressure 


D 


I 


I 


D 


D 


I 


I 


Mach 


I 


D 


D 


I 


I 


D 


to M<1 


Velocity 


I 


D 


D 


I 


I 


D 


D 


Temperature 


I/D 


I 


I/D 


D 


D 


I 


I 


Entropy 


I 


I 


D 


D 


I 


I 


I 


Stag. Press. 


D 


D 


I 


I 


D 


D 


D 


Stag. Temp. 


I 


I 


D 


D 








Density 


D 


I 


I 


D 


D 


I 


I 


Impulse Fct. 








1 D 


D 





Check this table with the information developed in the 
previous units and satisfy yourself that the Information 
is correct. 

h) Mow read Section 10.3 in the text. If time permits, 
glance at Sections 9-^, 9-5 and 10.5. 



175 



I 



i 



i 



1 

: 'I 



sr 

I 




i 




APPENDIX B 



STUDENT EVALUATION QUESTIONAIRE 

AE 30^3, "Gas Dynamics" - Quarter IV, 1971-72 



Please answer these questions as completely and candidly 
as possible. The results will in NO fashion effect your 
grade, but will be utilized in assessing the value and 
desirability of this method of instruction. 



V/here you are given a scale on vihich to indicate your opinion, 
a descriptive work (or v/ords) is provided to define each 
extreme. Please consider the scale linear with the mean 
midv;ay between 2 and 3- You may thus circle 1, 2, 3 or 4 
to indicate your own attitude relative to the extremes. 

Any amplifying comments you may wish to add to any question 
will be greatly appreciated. (Use the back of the question- 
aire or attach additional pages, as desired) 



1 . 



How many previous courses (of any duration, here at NPS 
or elsewhere) have you taken which were... 



SELF STUDY 
PROGRAMMED 



lAve. 1.9 



9 



SELF PACED 



75 P had had manyl ? 



lAve. 2.91 ? 



2. Indicate your opinion of the utility of the following 
educational innovations (NOTE: These features may or 



may NOT have 


been 


Included : 


in ' 


this 


course . ) 






a. Defined 


















Average 

Rating 


Obj ectives 


(no 


value ) 


1 


2 


3 




(great 


value ) 


3.50 


b. Self-Study 
Concept 


(no 


value ) 


1 


2 


3 




(great 


value ) 


1 

3.08 


c. Self-Paced 
Instruction 


(no 


value ) 


1 


2 


3 




(great 


value ) 


3.00 


d. Study Guide 


(no 


value ) 


1 


2 


3 




(great 


value ) 


3.67 1 


e. Check Tests 


(no 


value ) 


1 


2 


3 




(great 


value ) 


3.08 i 

1 


f. Group 

Discussions 


(no 


value ) 


1 


2 


3 




(great 


value ) 


1 

1 

3.18 


g. Individual 
Tutoring 


(no 


value ) 


1 


2 


3 




(great 


value ) 


1 

2.92 ' 


h. Audio- 
Tutorial 
(Tape-Slides ) 


(no 


value ) 


1 


2 


3 




(great 


value ) 


2.75 1 


i . Pinal Exam 


(no 


value ) 


1 


2 


3 




(great 


value ) 


2.67 



Boxed values indicate the average results of Questionnaires 
returned to the writer. 
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STUDENT EVALUATION QUESTIONAIRE (page 2.) 



The follov;ing questions pertain to the materials/methods 
used in this course. Please consider each question as it 
reflects on the course as a whole, then circle or underline 
your opinion. Then below each question, identify any unit 
(or section) to which your general answer does not apply and 
explain why. You v;ill probably find it useful to have your 
course materials at hand while considering these questions. 

3 . How would you rate the clarity of the Course OBJECTIVES? 

(obscure) 1 2 3 4 (clear) 

Exceptions : 



Average rating - 3-58 



Did you find the Study Guides understandable? 

(obscure) 123^ (clear) 

Exceptions : 

[Average rating - 3-51 



5 . What is your opinion of the organization of the Study 
Guides? 

(confused) 123^ (logical) 

Exceptions : 



Average rating - 3.^2 



6. Were sufficient graphs/diagrams/etc. provided in the 
Study Guides? 

(far too few) 123^ (far too many) 
Exceptions : 



Average rating - 2.^2 



7 . Were there enough example problems provided in the Study 
Guides? 

(far too few) 123^ (far too many) 

Examples : 



Average rating - 2.17 
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STUDENT EVALUATION QUESTIONAIRE (page 3-) 



8 . V/ere there enough homevjork problems provided to cover 
the material? 

(far too few) 123 ^ (far too many) 
Exceptions : 



Average rating - 3 -04 



9. V/as a correct portion of the quarter allocated for the 
completion of each unit In the course: 

(too little time) 1234 (too much time) 
Exceptions : 



Average rating - 2.58 



10. Outside of the scheduled problem sessions, on the average, 
how many times did you consult the Instructor for 
advice/assistance, per unit? 

Exceptions : 



Average - 1 time 



11. Did you seek advice/assistance from your classmates? 

(seldom) 1234 (regularly) 

Exceptions : 



Average rating - 2.5 



12. How do you feel the course materials aided you in pre- 
paring for the check tests? 

(useless) 1234 (Invaluable) 

Exceptions : 



Average rating - 3-58 



13 . How beneficial were the check tests in preparing for 
the major quizzes/final exam? 

(useless) 1234 (Invaluable) 

Exceptions : 



Average rating - 3-08 
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STUDENT EVALUATION QUESTIONAIRE (page H . ) 



1^. Estimate the percentage of your total study time 
expended on this course. 



Average rating - 25% 



15 . Do you feel any unit was excessively long or short? 
(If so please identify and comment.) 



Unit 


one - 


2 


Unit 


four- 


2 



16 . Hov; do you feel the amount you learned using this 

controlled, self-study method compares vjith v;hat you 
would have learned in a conventional lecture course? 

(less in this course; about the same; more in this course) 

0 

17* For several reasons this course V'jas run on a 
controlled -pace basis. 

a. Would you have liked the course better if it were 
completely self-paced? 



Yes 


- 5 


No 


- 6 



b. If it v/ere run on a completely self-paced basis, 

when do you think you would have finished the course? 



on time 


- 7 


early 


- 2 


late 


- 2 



18 . If you took a self-paced course and you were among the 
first students to successfully finish a unit, would 
you be willing to serve as a tuto r to those having 
more difficulty? (please commenFT 



Yes 


- 11 


No 


- 1 



19 . Assuming check tests would be retaken until mastery of 
material is essentially demonstrated, and that failure 
of a check test would not be held against you - v;ould 
it be better to have more comprehensive check tests and 
eliminate the major quizzes? (please comment) 



Yes 


- 6 


No 


- 6 
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STUDENT EVALUATION QUESTIONNAIRE (page 5.) 
20. V/ould you prefer check tests to be ... 



closed book? 

open book? 

part open, 
part closed 

oral? 



7 

0 

4 

1 



21. To what extent has participation in this course improved 

your ability to "decipher" a typical textbook or technical 
paper? 



little 8 

much 2 



22 . 



23. 



Given the choice of some type of self-study method as 
opposed to Instruction using the conventional lecture 
method, which v;ould you prefer? 

(You may qualify your answer in any way you choose... 
but please give your reasons . This is your chance 
to comment on anything you feel this questionnaire 
overlooked . ) 



self study 9 

lecture 2 



VJould you like to take all your courses by some type of 
a self-study method? 



no 

yes 

undecided 




THANK YOU FOR YOUR TIME AND EFFORT! 
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APPENDIX C 



SELECTED COMMENTS EXTRACTED FROM STUDENT QUESTIONNAIRES 
AND CRITIQUE SESSION. 



Question 19 

"... yes... eliminate the pressure of having to do 
v;ell at a specified time..." 

"Yes. Retakes under present check test, hour quiz policy 
are senseless waste of time." 

"No because it is very tiresome getting up for a test 
every week" 

"No. The check tests are aids, but the quizzes are for 
grades, (forgive my poetic license)" 



Question 21 



Comment by author; 

In general the students seemed to miss the intent 
of this question. Several commented that they had 
no basis on which to answer. Most answers seemed 
hesitant, or non-commital . 



Question 22 

"Self study Involving lesson guides, (audio-visual perhaps) 
Prob sets 

Followed by active participation seminar to reinforce 
important points" 

"Prefer a lecture. I have a good memory so I can keep 
up with previous ideas v;hile learning others. Most 
lecturers speak at an easily grasped pace for me." 



Comment by author: 

9 of the 12 questionnaires indicated a preference 
for this method or some modification of self-paced 
instruction . 



Comment from Critique Session 

"... This was really not v/hat I would call self-taught 
but was more like an ideal lecture course, where 
everyone has read his lesson assignment before class 
and the m.aterial is re-lnforced by the professor in the 
lecture ..." 
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